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A3STRACT 


Under  a three  year  program  investigations  were  conducted  on  the  solid 
state  reactions  which  occur  at  the  various  interfaces  in  thin  composite 
structures  on  the  elemental  and  compound  semiconductors,  silicon  and  gal- 
lium arsenide.  This  report  summarizes  the  results  of  these  investigations. 
Free  energies  derived  from  bulk  thermodynamic  data  are  used  to  predict  the 
direction  of  oxidation-reduction  reactions  at  interfaces.  The  identifica- 
tion of  the  resulting  solid  state  phases  is  accomplished  by  using  high 
energy  resolution  Auger  electron  spectroscopy  (AES)  and  X-ray  photoelectron 
spectroscopy  (XPS)  techniques  and  observing  the  chemical  shifts  in  energy 
associated  with  Auger  or  photoelectron  peaks  when  an  element  is  oxidized. 
Experimental  MOS  and  Schottky  barrier  test  structures  were  fabricated  in 
order  to  correlate  the  chemical  data  and  other  interface  characteristics 
with  electrical  measurements  made  on  these  devices.  Theoretical  interpre- 
tation of  the  results  of  high  energy  resolution  AES  and  XPS  studies  of  thin 
films  on  aluminum  gave  a value  for  the  attenuation  length  of  '1400  eV  elec- 
trons in  A1 203  and  information  on  various  electron  scattering  mechanisms. 

A review  of  GaAs  passivation  is  made  and  a method  based  upon  the  solid  state 
reactions  occurring  at  interfaces  is  suggested.  Finally,  localized  models 
for  predicting  surface/interface  states  on  the  elemental  and  compound  semi- 
conductors for  application  to  the  interface  problem  are  developed  and  photo- 
emission, energy  loss  spectroscopy  and  LEED  studies  are  examined  for  con- 
sistency with  the  predictions  of  the  theoretical  studies. 
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I.  INTRODUCTION: 


The  solid  state  device  requirements  for  advanced,  reliable  electronic 
systems  provide  a considerable  need  for  the  development  of  new  device  con- 
cepts and  fabrication  techniques  and  an  understanding  of  materials  inter- 
actions on  a fundamental  level.  The  purpose  of  the  research  program  de- 
scribed in  this  report  was  to  meet  this  need  in  part  by  investigating  the 
solid  state  reactions  which  occur  at  the  various  interfaces  in  device-type 
structures  on  the  elemental  and  compound  semiconductors,  silicon  and  gal- 
1 ium  arsenide. 

In  our  initial  proposal  of  February  18,  1975  we  presented  the  results 
of  calculations  of  the  standard  free  energy  change  on  the  formation  of  vari- 
ous inorganic  oxides  as  a function  of  temperature.  On  the  basis  of  these 
calculations,  that  is,  from  the  sign  of  the  free  energy  change  of  a solid 
state  reaction  between  two  solid  phases,  it  was  predicted  that  certain  metals 
on  silicon  dioxide,  for  example,  would  reduce  SiO^  and  form  a metal  oxide 
phase  at  the  interface.  In  order  to  identify  the  intermediate  oxide  phase 
it  was  proposed  to  use  the  techniques  of  Auger  electron  spectroscopy  (AES) 
and  X-ray  photoelectron  spectroscopy  (XPS)  taking  advantage  of  the  fact  that 
there  is  a chemical  shift  in  energy  between  the  peak  representing  the  metal 
and  that  representing  the  oxide.  During  the  first  year  preliminary  studies 
which  comprised  Auger  concentration  profiles  through  an  Al-SiO^-Si  structure 
and  XPS  data  at  the  Al/SiO^  interface  demonstrated  the  feasibility  of  iden- 
tifying the  intermediate  metal  oxide  phase  by  AES  and  XPS  techniques.  High 
energy  resolution  AES  and  XPS  studies  were  also  initiated  on  thin  oxide 
films  on  aluminum.  The  intensities  of  the  separate  peaks  for  A1  in  Al^ 
and  elemental  A1  were  determined  as  a function  of  oxide  thickness  and  the 
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direction  in  which  theoretical  interpretation  of  the  data  would  proceed  was 
outlined.  In  the  same  year  improvements  in  the  UHV/Auger  system  were  made 
and  an  evaporator  for  fabricating  MOS  or  Schottky  barrier  test  structures 
within  the  UHV  system  was  designed. 

During  the  second  year  the  theoretical  analysis  of  the  high  resolution 
Auger  and  XPS  data  for  the  thin  oxide  films  on  aluminum  was  completed  and 
an  attenuation  length  for  ' 1400  eV  electrons  in  AlpQ^  was  derived.  High 
resolution  Auger  and  XPS  data  for  oxidized  OaAs  surfaces  and  for  Al/GaAs  inter- 
faces were  also  obtained.  In  addition,  MOS  and  Schottky  barrier  planar  test 
structures  on  GaAs  were  prepared  in  a controlled  manner  in  the  Auger  system 
and  I-V  and  C-V  measurements  were  made  on  diodes  fabricated  from  the  planar 
structures  by  photoresist  and  masking  techniques. 

During  the  third  year  further  studies  were  pursued  on  the  preparation 
and  characterization  of  GaAs  surfaces,  particularly  on  molecular  beam  epi- 
taxial (MBE)  material,  and  additional  Schottky  barrier  devices  were  evaluated. 
The  problem  of  the  passivation  and  stabilization  of  GaAs  surfaces  was  then 
considered  as  an  application  of  the  fundamental  knowledge  of  solid  state 
reactions  at  interfaces.  The  final  effort  during  this  year  included  the 
beginning  of  a study  of  localized  theoretical  models  for  determining  surface/ 
interface  states  and  the  evaluation  of  existing  photoemission,  electron  energy 
loss  and  LEED  data  for  silicon  and  gallium  arsenide  surfaces. 

A discussion  of  the  thermodynamics  of  interface  reactions  is  presented 
in  Chapter  II.  The  basic  concepts  are  then  applied  in  studies  of  various 
interfaces  on  silicon  and  gallium  arsenide  (Chapters  III  and  IV).  The  high 
energy  resolution  AES  and  XPS  studies  of  Al^O^  on  A1  and  GaAs,  and  related 
optical  studies,  are  summarized  in  Chapters  V and  VI,  respectively.  GaAs 
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surface  passivation  is  covered  in  Chapter  VII  and  the  theoretical  study  of 
surface/ i nterface  states  is  detailed  in  Chapter  VIII.  Recommendations  are 
made  in  Chapter  XII  for  further  extensions  of  the  research  described  herein 
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I 1 • THERMODYNAMICS  OF  INTERFACE  REACTIONS: 

The  research  under  the  program  was  directed  towards  the  achievement 

of  a fundamental  understanding  of  the  physics  and  chemistry  of  interfaces 

on  elemental  and  compound  semiconductors.  One  important  problem  area  that 
I 

was  particularly  singled  out  was  the  nature  of  solid  state  reactions  oc- 
curring at  metal-oxide,  insulator-insulator,  oxide-semiconductor  and  metal- 
semiconductor  interfaces.  Thermodynamically  one  can  determine  the  likelihood 
of  the  chemical  reaction  that  might  occur  at  an  interface  between  two  solid 
phases  by  examining  the  standard  free  energy  change  of  the  process. 

The  Gibbs  free  energy  change  of  a reaction,  aG,  is  expressed  in  terms 
1 of  the  enthalpy,  AH,  and  the  entropy,  AS,  according  to  the  relation 

f 

AG  = AH  - TAS. 

The  coefficients  in  the  free  energy  equations  for  aG  and  aS  are  tabulated 

in  the  CRC  Handbook  for  Physics  and  Chemistry.  If  one  were  to  assume  that 

I 

aG  and  aS  were  independent  of  temperature,  then  the  temperature  dependence 
of  aG  could  be  obtained  readily  from  the  above  equation.  The  standard  free 
energy  change,  aG°,  refers  to  a reaction  for  which  all  reactants  and  pro- 
ducts are  in  their  standard  states  at  the  given  temperature  and  1 atm  pres- 
sure. Standard  free  energy  changes,  aG°,  for  the  formation  of  inorganic 
oxides  including  those  that  were  of  possible  interest  in  our  program  were 
' calculated  from  tabulated  data  for  each  comoound  for  a range  of  temperatures. 

The  results  are  shown  in  Figure  1. 

A graph  of  this  type  is  useful  in  evaluating  relative  reducing  or  oxi- 
dizing tendencies  of  materials.  For  example,  at  all  temperatures  A1  will 
reduce  Cr?0-i  to  Cr.  The  driving  force  for  this  oxidation-reduction  reaction 
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may  be  readily  obtained  from  consideration  of  aG°  for  the  two  compounds. 

At  600  °C 

2A1  (s ) + 3/2  02(g)  ■>  Al203(s)  aG°  = -355  kcal/mole 

2Cr(s)  + 3/2  0?  (g ) -»•  Cr203  (s)  aG°  = -234  kcal/mole 

Subtracting  the  second  reaction  from  the  first  yields 

2A1  (s)  +Cr203(s)  -*■  A1  ?03 (s ) + 2Cr(s)  aG°  = -121  kcal/mole 
The  negative  value  of  AG°  indicates  the  driving  force  for  the  reaction  re- 
ducing the  Cr203to  Cr  and  oxidizing  the  A1  to  form  A1203. 

Also  illustrated  in  Figure  1 are  differences  in  the  free  energy  change 
on  formation  for  V2 O5  and  Ga203,  at  700  °K,  and  for  A1203  and  Si02  at  800  °K. 
Consideration  of  the  aG°  for  the  first  pair  of  compounds  indicates  that  heat 
treatment  at  700  °K  of  vanadium  on  Ga203  will  result  in  the  formation  of 
V20ir  (or  other  vanadium  compounds  except  VO)  and  the  reduction  of  Ga203  to 
Ga  at  the  metal-oxide  interface.  Similar  arguments  apply  for  the  reduction 
of  Si02  by  aluminum  at  800  °K.  Note  also  that  according  to  the  bulk  data, 

Ti  is  expected  to  reduce  Ga203  at  high  temperatures  but  not  below  850  °K. 

A list  of  the  standard  heats  of  formation  (AH°)  and  standard  free  energies 
of  formation  (aG°)  of  various  oxides  is  given  in  Table  1.  Note  that  A1  is 
expected  to  reduce  Ga203  or  As203  to  form  A1203  and  free  Ga  or  As,  and  to 
reduce  Si 02  to  form  A1203  and  free  Si.  There  are  other  metals  which  are 
also  thermodynamical  ly  favorable  for  reducing  these  oxides,  e.g.,  Zr,  Cr, 

Nb  and  Ta. 

The  sign  of  the  free  energy  change  of  a solid  state  reaction  between 
two  solid  phases,  for  example,  the  metal  and  the  oxide,  of  course  governs 
only  the  feasibility  of  the  reaction.  Thus,  a aG  which  is  negative  will 
indicate  that  the  reaction  might  proceed.  The  kinetics  of  the  reaction. 
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Standard 

TABLE  I. 

Heat  of  Formation  (aH°)  and  Standard 

Free  Energy 

of 

Formation  (aG°)  of  Various  Oxides  at 

25°C 

(in  kilo-calories  per  gram  mole)* 

aH° 

aG° 

Ag202 

-6.62 

Ag20 

-7.74 

-2.59 

RhO 

-22.65 

Rh?0 

-23.74 

Cu20 

-39.84 

-34.98 

SeO^ 

-53.77 

Ru02 

-57.29 

CdO 

-60.86 

-53.79 

GeO 

-60.9 

SnO 

-68.4 

-61.5 

Rh20^ 

-70.06 

TeO 

-77.69 

-64.6 

Rb20 

-78.9 

Ga20 

-8T.0 

ZnO 

-83. T7 

-/6.05 

Rb202 

-101.7 

TcO^ 

-103.4 

TiO 

-125.0 

Ge02 

-128.3 

Tc03 

-129.0 

MoO? 

-130,0 

Sn02 

-138.8 

-124.2 

Cr03 

-141.0 

SrO 

-141.1 

-133.8 

Cr02  , 

-142.0 

MgO 

-143.84 

-136.13 

As203(Monoc1 inic) 

-150.76 

Sr02 

-153.6 

As20.,  (orthorhombic)  -154.87 

As?°4 

-173.69 
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TABLE  I. 


Standard  Heat  of  Formation  (aH°)  and  Standard  Free  Energy 
of  Formation  (aG°)  of  Various  Oxides  at  25°C 
(in  kilo-calories  per  gram  mole)* 

(Continued) 


AH° 

aG° 

MoO, 

yJ 

-180.33 

-161.95 

As2°5 

-218.6 

-184.6 

sb2o4 

-214.0 

-188 

Si02(quartz ) 

-205.4 

-192.4 

Mg (OH ) £ 

-221.0 

-199.27 

In2°3 

-222.5 

Ti02 

-228.0 

Sr(0H)2 

-229.3 

Sb2°5 

-234.4 

-200.5 

-258.0 

Zr02 

-258.2 

-244.4 

TcpOy 

-266.0 

Cr2°3 

-278.0 

Sb2°6 

-336.8 

-298.0 

y(oh)3 

-337.6 

-308.3 

Ti2°3 

-360.0 

NbpO^ 

-387.8 

Al203(corundum) 

-399.09 

-376.77 

Y2°3 

-419.6 

Nb205 

-463.2 

A 1 2O3 ’ ^2^ 

-471.0 

-435 

Ta2°5 

-493.0 

1 

•J. 

LO 

O 

-588.0 

Al203-3H20 

-613.7 

-547.9 

*from  Handbook  of  Chemistry  and  Physics, 
47th  Edition,  Chemical  Rubber  Co.,  (1966). 
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that  is,  the  reaction  rates  as  a function  of  temperature  need  to  be  deter- 
mined by  other  means.  There  are,  however,  other  factors  which  need  to  be 
considered  before  applying  these  thermodynamic  concepts  directly  to  the 
interface  reactions  on  the  elemental  and  compound  semiconductors.  The  first 
is  that  the  thin  film  situation  is  likely  to  be  somewhat  different  from  the 
bulk  so  that  the  free  energy  changes  might  not  be  accurate.  Secondly,  one 
should  be  concerned  about  the  distribution  of  the  elements  Ga,  As  or  Si 
after  the  reduction  reaction  has  occurred.  For  example,  it  is  expected 
that  these  elements  will  diffuse  from  the  interface  into  both  the  metal  and 

1 

1 2 

the  oxide  phases.  In  fact,  recent  Auger  studies  ’ on  the  interdiffusion 
between  aluminum  thin  films  and  GaAs  indicate  that  Ga  does  migrate  readily 
through  a 2000  A film  of  A1  at  250  °C. 

The  thermodynamic  concepts  presented  here  are  basic  to  the  studies  pre- 
sented in  later  chapters  and,  therefore,  reference  will  be  made  frequently 
to  this  chapter  and  the  expected  free  energy  change  of  a solid  state  reaction. 
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III.  INTERFACE  EFFECTS  ON  SCHOTTKY  BARRIERS: 


3 . 1 Introduction 

Several  experiments  were  performed  to  determine  the  effect  of 
interface  modifications  on  the  electrical  properties  of  Schottky  barriers 
formed  on  GaAs.  Surface  cleanliness,  stoichiometry,  lattice  damage  and 
oxidation  were  some  of  the  interface  conditions  investigated.  In  all  cases 
the  Schottky  barrier  structures  were  prepared  within  the  UHV/Auger  system. 
This  procedure  enabled  the  characterization  of  the  initial  GaAs  surface  and 
the  evaluation  of  cleaning  procedures  by  AES,  the  controlled  growth  of  an 
oxide  film  as  desired  with  monitoring  of  the  oxygen  (510  eV)  Auger  peak, 
and  deposition  of  metal  (Au  or  A1 ) using  a specially  designed  in  situ  eva- 
porator. Both  (100)  GaAs  substrate  material  and  (100)  GaAs  MBE  material 
were  used.  Subsequent  to  this  fabrication  procedure,  the  specimens  were 
removed  from  the  vacuum  system  and  actual  devices  were  formed  using  standard 
photoresist  and  masking  techniques.  Typically,  patterns  of  20  mil  dots  were 
prepared  on  the  metal  plating. 

The  evaluation  of  the  experimental  devices  involves  the  determination 

of  the  height  of  the  potential  barrier  formed  at  the  metal-semiconductor 

3 

interface  . There  are  several  techniques  that  can  be  used  to  determine  the 

4 

barrier  height  . The  instrumentation  we  used  enabled  measurements  of  the 
barrier  height  from  the  current-voltage  characteristics  and  the  capacitance- 
voltage  characteristics  of  the  Schottky  diodes. 

Current-voltage  characteristics  were  obtained  by  using  a Tektronics 
Type  575  transistor  curve  tracer  in  conjunction  with  an  external  picoammeter 
and  digital  voltmeter.  To  obtain  the  experimental  data,  the  current  through 
the  diode  was  measured  as  the  forward  bias  voltage  was  increased.  The  ideal 


10 


forward  I-V  characteristics  can  be  written  as 
J ■ A**T2  exp  exp 

where  <|>g  is  the  zero  field  barrier  height.  A**  is  the  effective  Richardson 
constant,  A<j>  is  the  Schottky  barrier  lowering,  V is  the  applied  voltage,  T 
is  the  temperature  and  k is  the  Boltzmann  constant.  If  the  saturation  cur- 
rent density  is  J$  at  zero  applied  bias,  then  the  barrier  height  can  be  ob- 
tained from 


By  plotting  J versus  V on  semi-log  paper,  the  saturation  current  may  be 
determined  by  extrapolation  of  the  current  density  at  zero  voltage. 

Capacitance-voltage  characteristics  were  measured  either  with  a Boonton 
Model  750  Direct  Capacitance  Bridge  or  the  Boonton  1 MHz  impedance  bridge 
which  is  an  integral  component  of  an  automatic  depth  profiler.  The  profiler 
measures  the  impurity  concentration  profiles  of  semiconductors  by  the  dif- 
ferential capacitance  method.  It  also  operates  in  the  C-V  mode,  whereby  it 
automatically  plots  the  capacitance  of  the  device  as  a function  of  applied 

4 

bias.  The  C-V  characteristics  are  given  by 

1/2 


where  q is  the  unit  charge,  is  the  doping  concentration,  e$  is  the  semi- 
conductor permitivity  and  the  built-in  voltage.  Thus,  from  the  voltage 

p 

intercept  of  a plot  of  1/C  versus  V,  the  barrier  height  can  be  determined 
from 
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Bn 


= V.  + V 
l r 


. kT  . , 


where  is  the  voltage  intercept  and  Vn  is  the  depth  of  the  Fermi-level 
below  the  conduction  band. 

3 . 2 Low  Pressure  Oxidation  of  GaAs 

The  formation  of  a natural  oxide  on  GaAs  under  low  oxygen  gas  pres- 
sure was  studied  in  one  series  of  experiments.  A small  wafer  of  (100) 
oriented  n-type  GaAs  substrate  material  was  obtained.  It  was  cut  into  two 

pieces  and  loaded  into  the  Auger  system.  The  system  was  pumped  overnight 

_g 

with  a 10  hour  bakeout  to  obtain  a pressure  of  1 x 10  Torr.  The  samples 
had  been  exposed  to  air  for  a long  time  and  therefore  initial  Auger  analysis 
showed  a great  deal  of  surface  carbon  and  oxygen.  One  specimen  was  then 
cleaned  by  argon  ion  sputtering  for  30  minutes  with  500  eV  ions  at  an  argon 

_5 

pressure  of  5 x 10  Torr.  The  specimen  was  then  examined  using  AES  and 

O 

then  plated  with  500  A of  pure  gold,  using  the  ijo  si tu  evaporator.  The 
Auger  spectrum  of  the  initial  surface  condition  before  Au  deposition  showed 
very  little  carbon  and  oxygen  remaining  on  the  surface  and  the  fine  structure 
in  the  low  energy  Ga  and  As  peaks  became  quite  apparent.  Since  the  low  and 
high  energy  peaks  of  Ga  and  As  have  different  escape  depths, s the  low  energy 
peaks  will  give  information  about  the  surface  composition  whereas  the  high 
energy  peaks  can  be  used  to  study  bulk  composition.  The  high  energy  peaks 

O 

of  Ga  (1070  eV)  and  As  (1228  eV)  have  escape  depths  of  19.0  and  20.5  A, 

O 

respectively.  The  escape  depth  is  5.0  A for  the  As  (31  eV)  peak.  No 
annealing  was  done  after  sputtering  so  that  lattice  damage  may  be  a contri- 
buting factor  to  the  results  obtained. 

The  second  specimen  was  cleaned  in  exactly  the  same  manner  and  showed 
almost  identical  surface  conditions  after  cleaning.  It  was  then  exposed  to 
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_5 

research  grade  oxygen  at  a pressure  of  1 x 10  Torr  by  backfilling  the  Auger 
system  using  the  automatic  pressure  controller.  The  low  energy  range  (20  to 
550  eV)  of  the  Auger  spectrum  was  monitored  as  a function  of  exposure  time. 
Particular  attention  was  paid  to  the  growth  of  the  oxygen  (510  eV)  peak. 

The  oxygen  peak-to-peak  height  grew  at  a constant  rate  of  .067  cm/min.  up 
to  120  minutes  exposure  time.  As  the  exposure  was  increased,  the  oxidation 
rate  tapered  off.  After  180  minutes  of  exposure,  the  oxygen  peak  had  in- 
creased from  .3  cm  to  12.8  cm.  At  this  point,  the  oxygen  was  pumped  out  and 

O 

500  A of  gold  was  deposited  on  the  sample. 

Afterwards,  both  specimens  were  removed  from  the  vacuum  system  and  20 
mil  diameter  dots  were  formed  on  the  plated  side.  Ni-Au  ohmic  contacts  were 
formed  on  the  reverse  side  of  each.  I-V  and  C-V  measurements  were  made  as 
described  earlier  and  barrier  heights  were  calculated.  From  the  I-V  measure- 
ments, barrier  heights  of  .85  volts  and  .93  volts  were  obtained  for  the  clean 
and  oxidized  specimen,  respectively.  C-V  measurements  gave  .84  and  .93  volts, 
respecti vely . Thus,  the  barrier  height  does  appear  to  depend  upon  the  inter- 
facial oxide  layer,  confirming  the  results  of  other  investigators,  Pruniaux 
c 

and  Adams,  for  example.  The  growth  rate  of  the  oxide  layer  in  our  case  was 

7 

probably  influenced  by  the  operating  hot  filament  ionization  gauge  as  well 
as  by  the  existence  of  lattice  damage  in  the  sputtered  GaAs  surface. 

3.3  Oxidation  of  GaAs  in  Air 

The  oxidation  of  GaAs  in  air  at  room  temperature  was  studied  and 
the  dependence  of  the  barrier  height  of  Schottky  diodes  with  interfacial 
layers  on  the  layer  thickness  was  determined.  High  resolution  Auger  and  XPS 
data  were  used  to  determine  the  nature  of  the  surface  oxide. 

A wafer  of  GaAs  was  obtained  and  cut  into  eight  pieces.  Ni-Au  ohmic 
contacts  were  formed  on  the  back  side  of  each  and  then  each  sample  was 
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chemically  cleaned  In  a sequence  of  methanol,  trichlorethylene  and  acetone 
baths.  Each  piece  was  then  etched  in  36%  HC1  for  20  seconds  and  kept  un- 
der methanol  until  loading  into  the  vacuum.  This  procedure  has  been  shown 

O 

to  leave  consistently  clean  surfaces  with  only  about  10  A of  oxide  on  the 

c _g 

surface.  The  system  was  pumped  down  overnight  to  a pressure  of  1 x 10 

Torr.  Initial  AES  examination  showed  the  presence  of  a great  deal  of  car- 
bon and  a small  Cl  peak. 

Three  of  the  specimens  were  ion  cleaned  to  give  surface  conditions 

O 

very  similar  to  those  described  in  Section  3.2.  About  500  A of  gold  was 

deposited  on  one  sample  as  is.  The  others  were  then  exposed  to  air  at  room 

temperature  for  one  and  two  hours,  respectively.  After  exposure  each  was 

O 

examined  by  AES  and  then  plated  with  500  A of  gold.  The  fourth  sample  was 

ion  cleaned  also,  but  much  more  thoroughly.  It  was  exposed  for  three  hours 

and  then  plated  in  the  same  manner. 

The  remaining  samples  were  ion  cleaned  to  give  conditions  similar  to 
the  first  three  specimens.  Two  were  oxidized  for  three  hours  in  air  and 
then  plated  with  A1  instead  of  Au.  The  remaining  sample  was  also  exposed 
for  three  hours,  but  no  metal  film  was  deposited  on  it.  All  the  samples 
except  for  one  of  the  A1  coated  specimens  and  the  unplated  specimen  went 
through  the  photoresist  process  to  form  diodes  for  the  barrier  height 
measurements . 

Figure  2 shows  a plot  of  the  barrier  height  versus  exposure  time  and 
the  oxygen  (510  eV)  peak-to-peak  amplitude  versus  exposure  time,  for  the 
first  four  samples.  Both  show  a linear  behavior  with  time  for  the  first 
three  specimens.  However,  the  fourth  sample  that  was  more  thoroughly 
cleaned  showed  a much  lower  oxidation  rate.  The  small  difference  in  the 
initial  surface  conditions  seems  to  have  quite  an  effect  on  the  oxidation 
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Figure  2.  Barrier  Height,  <t>gn,  of  Au/GaAs  Diodes  and  Oxygen  (510  eV) 
Peak-to-Peak  Amplitudes  versus  Exposure  Time  In  Air  at 
Room  Temperature  of  GaAs  Substrates  prior  to  Deposition  of 
Au  in  situ. 
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rate  and  the  barrier  height.  Pruniaux  and  Adams  also  found  two  different 
oxidation  rates  and  barrier  height  dependencies,  which  they  supposed  was 
caused  by  different  initial  surface  conditions,  although  they  had  no  way  of 
determining  the  surface  conditions  in  their  experiments. 

The  sample  which  was  oxidized,  but  not  plated,  was  used  to  obtain  high 
resolution  Auger  and  X-ray  photoelectron  spectra  (see  Section  5.2).  The 
oxide  on  the  GaAs  was  analyzed  before  and  after  sputter  removal  of  approxi- 

O 

rnately  15  A.  There  is  presently  some  disagreement  as  to  whether  the  oxygen 
on  the  surface  bonds  to  Ga  or  As  atoms  or  possibly  bridges  both.'"1'  Our 
high  resolution  data  indicate  the  presence  of  gallium  oxide  on  the  surface 
initially  and  no  arsenic  oxide.  See,  however,  the  results  of  photoemission 
studies  discussed  in  Chapter  VIII  which  indicate  that  oxygen  bonds  to  the  As 
atoms  on  a (110)  cleaved  surface. 

Barrier  height  measurements  were  made  on  the  Al/GaAs  diode  from  its 

? 

C-V  characteristics.  The  plot  of  1 /C  versus  V was  linear  over  the  voltage 
range  used.  Compared  to  the  reported  values  on  vacuum  cleaved  Al/GaAs 

4,11 

diodes  the  barrier  height  was  somewhat  larger.  The  high  resolution 
Auger  data  taken  on  the  other  sample,  which  was  prepared  in  exactly  the  same 
manner,  showed  no  evidence  of  either  gallium  oxide  or  an  arsenic  oxide. 

There  was,  however,  an  aluminum  oxide  at  the  interface,  extending  through  a 

O 

region  of  greater  than  55  A in  thickness  (see  Section  5.2).  These  results 
are  believed  to  be  direct  evidence  for  the  reduction  of  Ga^O^  by  A1  to  form 
Al^O^,  as  predicted  by  the  thermodynamic  considerations  discussed  in  the  pre- 
vious chapter. 

3 . 4 Experiments  with  MBE  Material 

A series  of  experiments  were  performed  using  (100)  GaAs  MBE  material 
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as  the  starting  material  for  Schottky  diodes.  A program  in  the  Solid  State 
Sciences  Division,  "Molecular  Beam  Epitaxial  Materials  Study  for  Microwave 
and  Millimeter  Wave  Devices,"  (Contract  No.  N001 73-76-C-0372) , provided  a 
ready  source  of  MBE  material  . 

following  the  MBE  run  the  GaAs  specimen  was  cut  in  half  and  both  pieces 

mounted  in  the  Auger  system,  one  on  a heated  stage.  The  system  was  then 

pumped  and  baked  for  approximatel y 16  hours  attaining  an  ultimate  pressure 
-9  ° 

of  2 x 10  Torr.  Approximately  500  A of  Au  was  deposited  on  one  specimen 
after  surface  characterization  by  AES.  The  other  specimen  was  cleaned  by 
argon  ion  sputtering  to  remove  the  carbon  followed  by  heating  to  600  °C  to 

O 

decompose  the  oxide.  Then  a 500  A film  was  deposited  ori  the  cleaned  sur- 
face. 

Low  energy  (10-650  eV)  and  high  energy  (650-1300  eV)  Auger  spectra  for 
the  specimens  showing  the  initial  surface  conditions  are  plotted  iri  Figure.. 

3A  and  3B,  respectively.  Note,  particular! y,  the  presence  of  carbon  at 
271  eV  and  oxygen  at  508  eV.  The  ratio  of  peak-to-peak  amplitudes  lor  the 
surface  sensitive  50  eV  (Oa)and  30  eV  (As)  peaks  is  1.6  while  the  ratio  for 
the  bulk  1066  eV  (Ga ) and  1227  eV  (As)  peaks  is  1.5. 

-5 

After  sputtering  in  argon  (0.5  keV,  25  ma,  5x1/1  Torr)  for  20  minutes 
the  low  and  high  energy  Auger  spectra  of  Figures  4A  and  4B  were  obtained. 

A considerable  reduction  in  the  carbon  peak  and  a marked  increase  in  the 
50  eV  (Oa)  and  30  eV  (As)  peaks  compared  to  the  same  peaks  in  Figure  3A  is 
noted.  The  ratio  for  the  50  eV  (Ga)  and  30  eV  (As)  peaks  is  0.8  while  that 
for  the  1066  eV  (Ga)  and  1227  eV  (As)  peaks  is  1.48,  approximately  the  same 
as  that  before  sputtering. 

Heating  to  600  °C  for  approx imately  4 minutes  resulted  in  the  complete 
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ELECTRON  ENERGY,  eV 


Figure  3.  Auger  Spectra  for  MBE  GaAs  Specimens 

with  Initial  Surface  Conditions.  A(10-650eV); 
B (650-1 300eV). 


Figure  4.  Auger  Spectra  for  MBE  GaAs  Specimens 

after  Argon  Ion  Sputtering.  A( 10-650eV) ; 
B (650-1 300eV). 
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disappearance  of  the  oxygen  peak  as  shown  in  Figure  5A,  however,  a depletion 

of  As  is  indicated  in  Figure  5B  by  the  change  in  the  ratio  of  the  peak-to- 

peak  amplitudes  for  the  1066  eV  (Ga)  and  the  1227  eV  (As)  peaks,  that  is, 

2.72  compared  to  the  previous  value  of  1.5. 

2 

Measurements  were  made  of  1/C  versus  V for  several  of  the  diodes  on 

the  wafer  which  had  gold  plated  over  the  existing  oxide.  Consistent  results 

were  obtained  showing  a linear  plot  and  an  intercept  in  the  neighborhood  of 
2 

1.1  eV.  The  1/C  vs  V curves  for  the  cleaned  GaAs,  however,  showed  a non- 
linearity at  higher  voltages  and  a range  of  intercept  values  from  0.0  to 

1.2  eV.  This  nonreproducibility  is  possibly  associated  with  the  nonuniform- 
ity of  cleaning  and  the  nonlinearity  with  the  As  depletion.  Further  studies 
of  this  type  are  indicated  to  define  more  clearly  the  relation  between  the 
particular  surface  conditions  and  the  capacitance-voltage  characteristics. 
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Figure  5.  Auger  Spectra  for  MBE  GaAs  Specimens 
after  Heating  to  600°C.  A(10-650eV); 
B(650-l 300eV ) . 
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IV.  INTERFACE  STUDIES  OF  A1 -Si02~Si  STRUCTURES: 

In  preliminary  work  on  MOS  structures  fabricated  by  depositing  A1  on 
thermally  oxidized  silicon  wafers,  as  reported  in  our  Interim  Technical 
Report  of  21  January  1976,  it  was  shown  that  the  chemical  shift  in  energy 
between  the  Al-2p  photoelectron  peak  due  to  A1  and  that  due  to  A1  in  AlpO^ 
was  sufficient  to  identify  the  two  separate  solid  state  phases  as  one  sput- 
tered through  the  Al/SiO^  interface  of  an  Al-SiO^-Si  structure.  The  spectra 
obtained  indicated  that  an  Al^O^  layer  was  formed  at  the  Al/SiO^  interface 
according  to  the  predictions  of  bulk  thermodynamics , as  discussed  in  Chapter 
II.  The  Si  peaks  were  not  monitored  at  that  time  so  that  the  results  were 
inconclusive  with  respect  to  the  reduction  of  SiO^  by  the  A1  to  form  free 
silicon. 

Subsequently,  similar  specimens  were  examined  more  thoroughly  in  a 
Varian  Auger  electron  spectrometer  which  incorporated  a cylindrical  mirror 
analyzer  with  a 0.25%  energy  resolution.  Auger  derivative  data  were  taken 
using  four  channels  of  a multiplexer  to  obtain  an  Auger  profile  thorugh  the 
Al-SiO^-Si  structure.  The  lower  channel  included  both  the  Al-LVV  peaks  due 
to  the  A1  and  Al^O^  and  the  Si-LVV  peaks  due  to  Si  and  SiO^;  the  second 
channel  included  only  the  0-KLL  peak;  the  third  channel  included  the  KLL 
peaks  due  to  A1  and  Al^O^;  and  the  fourth  channel  included  the  KLL  peaks  due 
to  Si  and  SiO,,.  Since  the  Auger  profile  is  a point  plot  of  the  largest  peak 
in  each  channel  as  a function  of  sputtering  time  it  is  not  possible  to  dis- 
tinguish between  A1  and  AlpO^  or  Si  and  SiO^  peaks  from  such  a profile. 
However,  we  have  simultaneously  plotted  the  data  within  each  channel  on  a 

strip  chart  recorder.  From  the  strip  chart  records  one  can,  therefore,  ob- 
tain information  on  the  concentrations  of  the  different  phases  through  the 
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Al/SiOp  interface. 

The  striking  results  of  these  studies  show  the  presence  of  free  A1 , 

A1203,  free  Si  and  SiO^  at  the  Al/SiO^  interface;  an  A^O^  distribution 

o 

extending  from  the  interface  approximately  150  A into  the  A1  and  approxima- 

o 

tely  150  A into  the  Si02;  and  a distribution  of  free  Si  extending  partly 

o o 

into  the  A1  (<  150  A)  and  to  a greater  extent  into  the  Si02  (>  200  A).  The 
broad  reaction  zone  was  not  expected  in  that  the  sample  was  maintained  at  a 
relatively  low  temperature  (:;  70  °C  or  less)  during  the  A1  deposition.  How- 
ever, the  energy  released  in  the  solid  state  reaction  involving  the  reduction 

1 

of  Si02  and  the  formation  of  A1203  could  be  sufficient  to  cause  considerable 
intermixing  of  the  several  components.  A consequence  of  these  studies  is 
that  further  annealing  of  the  A1  films  as  is  ordinarily  done  in  the  produc- 
tion of  M0S  devices  (typically,  470  °C  in  dry  N2  for  20  min.)  will  probably  not  ; 

appreciably  extend  the  chemical  reduction  of  Si02  by  A1 . Further  studies 
of  the  phenomena  involving  different  deposition  methods  and  conditions, 
different  metals  and  thinner  oxide  films  are  indicated. 

A paper  covering  the  details  of  these  experiments  by  Y.  F.  Strausser 
(Varian),  E.  J.  Scheibner  and  J.  S.  Johannessen  (Stanford)  has  been  accepted 
for  publication  in  Thin  Sol  id  Films. 
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V.  HIGH  RESOLUTION  AUGER  AND  X-RAY  PHOTOELECTRON  SPECTROSCOPY  OF  ALUMINUM 
AND  GALLIUM  ARSENIDE  SURFACES: 

In  this  phase  of  the  research  a double-pass  cylindrical  mirror  analyzer 

] 2 

with  retarding  grids  was  used  to  obtain  high  energy  resolution  electron- 
excited  Auger  and  X-ray  photoelectron  spectra  from  thin  oxide  films  on 
aluminum  and  from  different  oxidized  gallium  arsenide  surfaces.  With  the 
double-pass  CMA  the  incident  Auger  electrons  (or  photoelectrons)  are  de- 
celerated in  the  retarding  field  of  spherically-shaped  grids  in  such  a manner 
that  a constant  pass  energy  through  the  CMA  is  obtained.  The  result  is  that 
the  absolute  resolution,  as  indicated  by  the  FWHM  of  the  elastic  peak,  is 
a function  only  of  the  pass  energy  and  the  size  apertures  used  in  the  analy- 
zer and  is  not  dependent  upon  the  incident  electron  energy,  as  it  is  for 
the  single-pass  CMA.  For  the  conditions  we  used,  that  is,  pass  energies  of 
the  order  of  50  eV  and  small  apertures,  the  absolute  resolution  was  about 
0.7  eV.  This  made  it  possible,  for  example,  to  distinguish  easily  between 
the  Auger  KLL  peaks  due  to  A1  in  Al^O^  and  those  due  to  A1  metal,  which  are 
separated  in  energy  by  about  6.7  eV.  In  fact,  both  Auger  and  X-ray  photo- 
electron spectra  exhibit  such  chemical  shifts  for  guite  a number  of  metal - 
1 1 

metal  oxide  pairs.  In  general,  the  Auger  chemical  shifts  are  considerably 
larger  than  the  photoelectron  shifts  when  the  transitions  involve  only  core- 
level electrons. 

The  analysis  of  the  Auger  KLL  spectra  and  the  Al-2p  photoelectron  spectra 
for  different  o/.ide  film  thicknesses  on  aluminum  films  was  presented  in  a 
paper  submitted  to  the  Journal  of  Electron  Spectroscopy  and  Related  Phenomena. 
A copy  of  that  paper  is  also  included  in  the  Interim  Technical  Report  dated 

1 4 

28  January  1977.  Studies  of  the  gallium  arsenide  surfaces  are  also  included 
in  the  same  report.  The  principal  results  of  both  studies  will  be  summarized 
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here. 

5 . 1 Aluminum/Aluminum  Oxide 

Specimens  used  in  this  investigation  were  aluminum  films  sputtered 
in  a Sloan  DC  Sputtergun  system  onto  thermally  oxidized  silicon  wafers.  The 
aluminum  oxide  films  were  the  natural  oxide  of  aluminum  formed  in  air  during 
24  hours  in  the  sputtering  system.  Sputter  etching  of  the  Al^O^  films  to 
obtain  different  oxide  film  thicknesses  was  accomplished  by  argon  ion  bom- 

O 

bardment  in  the  Auger/XPS  system  at  a rate  of  approximately  30  A/min.  Auger 
and  XPS  spectra  were  recorded  initially  and  after  each  stage  of  argon  ion 
sputtering. 

Figure  6 is  a high  resolution  Auger  spectrum  obtained  from  the  surface 
of  an  oxidized  aluminum  specimen  showing  the  major  peaks  in  the  A1  KLL  spec- 
trum. These  include  the  most  intense  KLL  peaks,  the  KLpL-^D^)  and  the 
KL-^^P-j)  peaks,  and  the  weaker  KLgLgx  Sg)  and  KL-jL^^Pq  1 ,,)  peaks.  Addi- 
tional KLL  peaks  should  occur  at  approximately  1301  eV  for  the  KL^L^(^Sq) 
transition,  which  is  outside  the  range  of  our  data,  and  at  1392  eV  for  the 
KL^L3(  Pq  2)  transitions,  which,  however,  is  below  the  limit  of  detection 
of  our  instrument.  Bulk  plasmon  loss  peaks  associated  with  the  KL^L^D) 
and  the  KL-^^P)  Auger  electrons  as  well  as  surface  plasmon  peaks  associated 
with  the  KL^L^'d)  Auger  electrons  are  indicated.  The  lower  curves  in  the 
figure  which  are  energy  loss  data  for  a primary  electron  energy  correspond- 
ing to  the  KL^L-^D,,)  Auger  energy  show  the  bulk  plasmon  loss  peaks  in  mul- 
tiples of  15.6  eV  from  the  primary  energy.  The  surface  plasmon  loss  peaks 
in  the  Auger  spectrum  occur  at  a Ep//T  = 11.0  eV  and  at  E //T  = 9.0  eV 
from  the  KL^L^D^)  Auger  energy.  These  values  are  related  to  the  metal- 
vacuum  interface  topography,  a planar  interface  giving  an  energy  difference, 


25 


KLj-.fD,) 


ELECTRON  ENERGY.  eV 


Figure  6.  High  Resolution  Auger  Spectrum  from  Aluminum  Showing 

KLL  Transitions  and  Bulk  and  Surface  Plasmon  Loss  Peaks 
The  lower  curves  are  electron  energy  loss  data  for 
primary  electrons  of  the  KL2L3(’D2)  Auger  energy. 


aE,  equal  to  Ep//T  and  a rough  or  spherical  interface  giving  a aE  of  E //T. 
Figure  7 shows  high  resolution  Auger  spectra,  N(E)  vs  E,  for  the  KL^ 

o 

transitions  before  argon  ion  sputtering  (curve  A),  after  15  A removal  (curve 

o o 

B),  after  30  A removal  (curve  C)  and  for  the  aluminum  at  75  A removal  (curve 

o 

D).  The  total  thickness  of  the  oxide  film  was  of  the  order  of  40  A.  Large 


Auger  chemical  shifts  of  over  7 eV  between  the  metal  and  oxide  peaks  are 
shown  in  curves  A and  B.  The  slight  differences  may  be  due  to  compositional 
variations  in  the  oxide  or  to  surface  charging.  The  aluminum  curve  (curve 
D)  shows  weak  peaks  at  1378.3  eV  and  at  1382.1  eV.  By  comparison  with  Fi- 
gure 6 the  1382.1  eV  peak  is  identified  as  the  KL^L^Sq)  Auger  peak  and 
the  1378.8  eV  peak  is  associated  with  the  excitation  of  surface  plasmons  by 
Auger  electrons  from  the  K^L^D^)  peak.  The  increasing  intensity  at  the 
left  of  each  curve  results  from  the  bulk  plasmon  loss  peak.  Curve  C includes 
contributions  from  the  KL^L^Sq)  Auger  peak  and  from  the  surface  plasmon 
peak.  In  our  analysis  we  find  that  the  surface  plasmon  peak  shifts  to  a 
lower  aE  value  in  the  presence  of  a thin  oxide  layer  as  suggested  in  the  fi- 
gure. For  a thick  oxide  the  appropriate  relation  for  the  surface  plasmon 

1 5 

energy  would  be  E //  1+e  where  e is  the  dielectric  constant  of  the  oxide. 

High  resolution  photoelectron  spectra  of  the  Al-2p  photoelectron  peak 
are  shown  in  Figure  8 for  an  oxidized  aluminum  specimen  before  sputtering 

O 

(curve  A),  after  30  A removal  (curve  B)  and  for  the  aluminum  (curve  C). 

The  smaller  chemical  shift  between  the  Al-2p  aluminum  peak  at  71.7  eV  of 

1 3 

2.5  eV  is  in  agreement  with  Wagner's  data 

The  different  Auger  and  X-ray  photoelectron  spectra  were  analyzed  as 
follows.  The  experimental  shape  of  a single  peak  is  usually  a Voigt  func- 
tion, a convolution  of  Lorentzian  (L)  and  Gaussian  (G)  functions.  Ue 
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Figure  7.  Auger  Spectra  for  the  KL2L3  Transition  Showing  Metal  and  Oxide 
Peaks  for  Different  Thicknesses  of  A1 .0?  on  Aluminum.  Oxide 
thickness  is  of  the  order  of  40  A.  Curve  A is  for  the  specimen 
before  argon  ion  sputtering,  curve  B is  after  15  A removal, 
curve  C is  after  30  A removal  and  curve  D is  for  the  aluminum 
without  an  oxide. 


?8 
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approximated  the  Voigt  function  by  a linear  combination  of  L and  G compo- 
nents involving  the  height,  position  and  FWHM  of  the  peak  and  a L-G  mixing 

i c 

coefficient.  We  then  fitted  all  the  peaks,  one  peak  at  a time,  using  a 
non-linear  least  squares  fitting  procedure.  A background  subtraction  which 
includes  the  characteristic  loss  spectrum  is  made  and  an  iteration  procedure 
is  followed  until  the  best  fit  is  obtained.  The  results  of  the  fitting  pro- 
cedure for  the  Auger  data  are  shown  in  Figure  9 where  the  solid  lines  are 
the  experimental  data.  For  the  XPS  data,  it  should  be  noted  that  each  peak 
is  an  unresolved  doublet  involving  the  2p-,^  and  Pp^,,  states.  With  the 

1 7 

theoretical  ?p.^/2p-|/p  intensity  ratio  of  2 and  a splitting  of  0.4  eV, 
the  peaks  in  curves  A and  B were  successfully  fitted.  A "relaxed"  fit  was 
obtained  for  curve  C by  allowing  the  2p  doublet  to  relax  for  the  best  fit 
without  imposing  restrictions  on  the  splitting  or  the  intensity  ratio.  The 
results  of  the  fitting  procedure  for  the  XPS  data  are  shown  in  Fiqure  10. 

After  the  fitting  was  completed  for  both  Auger  and  XPS  peaks  and  the 
contribution  from  the  background  and  neighboring  peaks  removed,  the  area 
under  each  metal  and  oxide  peak  was  found  by  integrating  under  the  experi- 
mental curve.  The  attenuation  length  x,  in  AlpO^,  for  electrons  of  approxi- 
mately 1400  eV  was  then  found  by  computing  the  effect  of  geometry  and  the 
oxide  layer  thickness  on  the  Auger  and  XPS  intensities  and  including  appro- 
priate emission  functions  for  the  oxide  and  the  metal  as  described  by 

l a o 

Shelton.  An  average  value  of  x = 31  + 3 A was  obtained  by  taking  the  mean 
of  all  the  attenuations  for  the  metal  and  oxide  peaks  in  both  the  AES  and 
XPS  spectra.  The  oxide  and  metal  normalized  integrated  intensities  plotted 
versus  t/x  are  shown  in  Figure  11. 


30 


N(E). ARBITRARY  UNITS 


BINDING  ENERGY.  eV 

Figure  10.  Computer  Results  of  Fitting  Theoretical  Functions  to  Experimental  XPS  Spectra. 

Solid  lines  are  smoothed  experimental  data;  dots  are  computed  points.  The 
designated  curves  correspond  to  those  of  Figure  8. 
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Attenuation  Plots  Showing  Nornallzod  Integrated  Intensities  versus  Oxide  Thickness  for 
the  Metal  and  Oxide  Peaks  from  AES  and  XPS  Spectra. 


The  value  of  the  attenuation  length  may  be  in  error  due  to  the  difference 
in  sputtering  rate  for  the  calibration  film  used,  Ta^O^,  and  Al^O^.  Further 
experiments  during  the  course  of  the  program  were  directed  toward  the  pre- 
paration of  thin  oxide  films  on  tantalum  and  aluminum  by  anodization  and 
reactive  sputtering,  and  an  independent  measurement  of  the  film  thickness 
by  el  1 ipsometry.  Using  such  specimens  with  known  film  thicknesses  in  each 
Auger  or  XP5  run  would  enable  a more  precise  determination  of  the  sputtering 
rates.  When  these  experiments  are  completed  a reevaluation  of  the  attenua- 
tion data  is  anticipated.  A summary  discussion  of  the  optical  methods  used 
in  the  program  is  included  in  the  next  chapter  of  this  report. 

The  results  of  the  computer  analysis  also  yield  information  on  lifetime 
and  phonon  broadening.  These  topics  are  discussed  fully  in  the  complete 
paper 

5.2  Gallium  Arsenide 

Three  types  of  GaAs  samples  were  prepared  and  examined  by  high 
resolution  Auger  and  XPS  techniques:  (1)  a (100)  GaAs  specimen  oxidized  in 
air  at  room  temperature;  (2)  an  Al/GaAs  junction  formed  by  depositing  A1 
in  situ  on  a (100)  GaAs  specimen  oxidized  in  air  for  3 hours  at  room  tempera- 
ture; and  (3)  a (100)  GaAs  specimen  thermally  oxidized  in  0^  at  600  °C  for 
20  minutes. 

Figure  12  shows  the  Ga  2p ^ photoelectron  peaks  for  a gallium  arsenide 
sample  oxidized  in  air  at  room  temperature.  Curve  A is  for  the  oxidized 

o 

specimen;  curve  B is  for  the  same  specimen  after  removing  approximately  15  A 
by  argon  ion  sputtering.  A chemical  shift  in  energy  between  the  peak  due  to 
gallium  and  the  peak  due  to  gallium  in  the  oxide  is  about  1.2  eV.  It  should 
be  noted  also  that  because  of  the  high  binding  energies  (1117  and  1118.2  eV) 
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Figure  12.  Ga  2pi/?  Photoelectron  Peaks  for  Gallium  Arsenide  Oxidized 
in  Air  at  Room  Temperature.  Curve  A,  showing  a chemical 
shift  of  about  1.2  eV,  is  for  the  oxidized  specimen;  curve 
B Is  for  the  same  specimen  after  removing  approximately  15  A 
by  argon  ion  sputtering.  The  number  of  scans  in  the  signal 
averager  is  indicated  on  each  curve. 


and  the  MgKa  photon  energy  of  1253.6  eV  the  kinetic  energies  of  the  photo- 
electrons are  only  of  the  order  of  130  eV.  The  escape  depth  for  electrons 
of  this  energy  is  small  enough  that  there  is  practically  no  contribution 
to  the  oxide  peak  due  to  photoelectrons  from  the  Ga  in  the  substrate. 

Figure  13  shows  the  X-ray  excited  Auger  peaks  for  the  same  gallium 
arsenide  specimen  which  was  oxidized  in  air  at  room  temperature.  Curve  A 
is  for  the  oxidized  specimen;  curve  B is  for  the  same  specimen  after  sputter 

O 

removal  of  15  A.  The  Auger  peaks  occur  in  the  XPS  spectrum  at  apparent 
binding  energies  equal  to  the  photon  energy  less  the  sum  of  the  Auger  energy 
and  the  spectrometer  work  function.  For  these  peaks  the  Auger  energies  are 
about  1060  eV  and  there  will  be  an  appreciable  contribution  to  the  Auger 
peaks  due  to  the  oxide  from  the  Auger  electrons  of  the  Ga  in  the  substrate 
because  of  the  greater  escape  depth  at  the  higher  energy.  The  spectra  in 
Figure  12  and  13  show  that  room  temperature  oxidation  produces  a gallium 
oxide;  no  evidence  for  an  arsenic  oxide  was  found  in  either  Auger  or  XPS 
data  for  the  specimens  examined.  However,  a different  pre-oxidation  surface 
treatment  and/or  a different  crystallographic  orientation  may  show  the  growth 
of  arsenic  oxide  at  room  temperature.  Detailed  oxidation  studies  would  be 
required  to  establish  the  conditions  for  the  growth  of  either  or  both  type 
oxide  films. 

In  order  to  determine  the  effect  of  an  A1  layer  deposited  on  the  room 
temperature  oxide  a GaAs  specimen  was  oxidized  in  air  for  3 hours  at  room 
temperature  and  then  an  A1  layer  was  deposited  in  the  ultra-high  vacuum 
Auger  system  using  an  in  situ  evaporator.  A similar  specimen  was  prepared 
for  the  Schottky  diode  measurements  reported  in  chapter  III.  Figure  14 
shows  the  Al-KLL  Auger  spectra  for  different  thicknesses  through  the  Al/GaAs 
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Figure  13.  X-ray  Excited  Ca  Auger  Peaks  for  Gallium  Arsenide  Oxidized 
in  Air  at  Room  Temperature.  Curve  A is  for  the  oxidized 
specimen,  curvg  B is  for  the  same  specimen  after  sputter 
removal  of  15  A.  The  number  of  scans  is  indicated  on  each 
curve.  Auger  peaks  occur  in  the  XPS  spectrum  at  apparent 
binding  energies  equal  to  the  photon  energy  (for  MgKa  radiation, 
the  photon  energy  is  1253.6  eV)  less  the  sum  of  the  Auger 
energy  and  the  spectrometer  work  function.  For  example, 

B.E . = 1253.6  - (1058  ♦ 4)  = 191.6  eV. 
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Mf’ure  14.  A7.-KLL  Auger  Spectra  for  Different  Thicknesses  through  the 
Interface  of  an  At/GaAs  Junction  formed  by  Depositing  Af 
in  situ  on  a GaAs  Specimen  Oxidized  in  Air  for  3 Hours  at 
Room  Temperature.  Curve  A is  after  removal  of  approximately 
130  A of  A?.;  curve  B is  after  10  A removed  from  curve  A depth; 
curve  D is  after  40  A removed  from  curve  A depth;  and  curve  F 
is  after  55  A removed  from  curve  A depth.  The  number  of  scans 
is  indicated  on  each  curve. 


interface.  Curve  A is  after  removal  of  approximately  130  A of  A1  ; curve  B 

O O 

is  after  10  A removed  from  the  curve  A depth;  curve  C is  after  20  A removed 

O 

from  the  curve  A depth;  curve  D is  after  40  A removed  from  the  curve  A 
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depth;  and  curve  E is  after  55  A removed  from  the  curve  A depth.  The  pre- 
sence of  A1  and  A1,,0,,  in  a broad  interface  region  is  similar  to  the  result 
obtained  for  the  Al/SiO^  interface  (Chapter  IV).  Arsenic  and  gallium  LMM 
Auger  spectra  showed  only  elemental  As  and  Ga,  and  no  residual  Ga^O^  or  As^O^. 

The  third  sample  examined  was  thermally  oxidized  GaAs  specimen.  The 
purpose  was  to  determine  the  composition  of  the  oxide,  mainly  near  the  outer 
oxide  surface.  Figure  15  shows  Ga  2p.^?  photoelectron  peaks  for  the  oxide 
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surface  (curve  A);  after  30  A removed  (curve  B);  and  after  60  A removed 
(curve  C).  The  interpretation  is  that  the  oxide  is  Ga^O^.  The  peak  posi- 
tion agrees  with  the  oxide  peak  position  in  Figure  12.  Note  that  the  peak 
intensity  increases  then  remains  relatively  constant  through  the  oxide  film. 
Figure  16  shows  As  3d  photoelectron  peaks  for  the  oxide  surface  (curve 
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A);  after  30  A removed  (curve  B);  after  60  A removed  (curve  C);  and  after 
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120  A removed  (curve  D).  The  As  peak  at  41.5  eV  builds  up  to  a constant 
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value  after  removing  approximately  120  A while  the  peak  due  to  As  in  As^O^ 
decreases  rapidly  from  a high  value  at  the  surface.  The  interpretation  is 
that  there  is  an  As^O^  7ayer  at  the  surface  which  decreases  exponentially 
into  the  Ga^O^,  and  that  there  is  free  As  distributed  through  the  Ga^O^ 
thickness.  The  lower  Ga^O^  peak  intensity  for  curve  A of  Figure  15  is  con- 
sistent with  this  picture  since  the  Auger  electrons  from  the  Ga  in  Ga^O^ 
would  be  attenuated  by  the  As^O^  layer. 

The  significant  result  of  this  section  is  the  apparent  reduction  of 
Ga^Oj  by  A1  to  form  Al^O^  at  the  interface.  The  stability  of  Schottky  diodes 
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involving  other  reactive  metals  as  well  as  aluminum  and  the  possibility  of 
utilizing  such  solid  state  reactions  to  develop  a passivation  method  for 
GaAs  devices  have  been  considered  in  our  program.  Previous  passivation 
methods  are  reviewed  in  chapter  VII  followed  by  ari  outline  of  a suggested 
method  based  upon  the  reduction  of  the  initial  gallium  or  arsenic  oxide. 
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VI.  OPTICAL  METHODS: 


In  parallel  with  the  previous  studies  ari  effort  was  initiated  to  charac- 
terize thin  oxide  films  on  gallium  arsenide  and  other  materials  by  the  op- 
tical methods  of  ellipsometry  and  reflectivity.  Oxide  films  on  aluminum 
and  tantalum,  prepared  by  anodic  oxidation  or  reactive  sputtering,  were  studied 
to  determine  the  optical  constants  and  film  thickness.  With  a series  of 
these  films  with  calibrated  thicknesses  one  could  in  principle  determine 
sputtering  rates  and  thus  obtain  a more  accurate  value  for  the  electron  at- 
tenuation lengths  in  oxides.  In  fact,  it  would  be  desirable  to  incorporate 
a calibration  specimen  along  with  the  specimen  being  examined  by  high  reso- 
lution AES  and  XPS  in  order  to  insure  that  the  sputtering  conditions  were 
the  same  for  both  type  specimens. 

The  ellipsometer  used  in  our  investigations  was  a LI  17  Gaertner  Pro- 

o 

duction  Ellipsometer  equipped  with  a 2 mW  helium-neon  laser  (6323A)  light 
source.  The  el  1 ipsometric  parameters  that  are  measured  are  A and  ip . A is 
defined  as  6 -■!'c,  the  difference  between  the  phase  shifts  for  the  component 

r ^ 

of  the  light  wave  parallel  (p)  to  the  plane  of  incidence  and  that  perpendi- 
cular (s)  to  the  plane  of  incidence.  The  relative  amplitude  attenuation  is 
given  by  tan  >p.  A and  are  related  to  the  optical  constants  through  the 
relation 

tan  <|>  e1A  = r /r 

r 

where  rp  and  r^  are  the  Fresnel  reflection  coefficients  for  a surface  film- 
absorbing substrate  system.  For  such  a system  one  obtains  an  expression 
relating  five  known  parameters  (> , <{> , n^ , A and  to  five  unknown  parameters 
(n^,  kp,  d,  n^,  and  k^).  We  have  used  the  computer  program  of  Mai  in  and 
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Vidam  to  determine  the  optical  constants  for  the  thin  film  and  absorbing 

substrate  and  the  thickness  of  the  oxide  film  from  four  sets  of  A and  <p 

? c 

values.  We  have  also  used  the  computer  program  of  So  t.o  determine  the 
optical  constants  of  the  film  and  its  thickness  from  only  two  sets  of  A 
and  ip  values  when  the  optical  constants  of  the  substrate  are  known.  Ob- 
taining the  optical  constants  of  a clean  substrate  is  not  a trivial  matter. 

In  our  studies  we  used  for  this  purpose  an  automated  surface  analysis  appa- 

.■  1 

ratus  developed  by  Wall  and  Stevenson  of  Georgia  Tech  which  provides  real 
time  acquisition  and  manipulation  of  optical  reflectivity  and  Auger  data. 
The  apparatus  consists  of  an  ultra-high  vacuum  chamber/ ref lectometer  equip- 
ped with  a cylindrical  mirror  analyzer  and  integral  coaxial  gun  for  surface 
characterization  by  AES,  an  ion  gun  tor  surface  cleaning  by  argon  ion  sput- 
tering and  a minicomputer  for  the  collection  and  manipulation  of  the  data. 
Computer  programs  were  available  to  determine,  from  the  reflectivity  data, 
the  optical  constants,  n arid  k,  for  the  range  of  wavelengths  used. 

Gallium  arsenide  substrates  were  examined  by  both  optical  methods. 

The  thickness  of  thin  oxide  films  on  GaAs  after  wet  chemical  cleaning  was 
determined  by  ellipsoinetry  and  reflectivity  measurements  for  clean  and  ano- 
dized GaAs  were  compared. 


i 
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VII.  GaAs  SURFACE  PASSIVATION  - A PROPOSED  METHOD: 


Significant  applications  of  GaAs  devices  and  integrated  circuits  depend 
upon  the  development  of  an  oxide  or  insulating  layer  v/i th  good  dielectric  pro- 
perties, a minimum  of  fixed  and  mobile  ions  iri  the  layer,  and  low  surface  or 
interface  state  densities  at  the  insulator-GaAs  interface.  In  the  course  of 
our  fundamental  investigations  on  GaAs  surfaces  we  have  also  been  cognizant 
of  this  need  for  surface  passivation  and  have  considered  the  application  of 
our  results  to  the  passivation  problem.  A brief  survey  of  the  literature 
was  made  and  a practical  method  was  then  outlined  which  utilizes  the  principle 
of  the  solid  state  reduction  of  an  initial  oxide  on  GaAs  by  A1  to  form  AlpfJ^. 

A more  detailed  study  would  of  course  be  necessary  to  explore  this  approach 
but  the  understanding  of  the  surface  physics  and  chemistry  at  insulator-GaAs 
interfaces  is  believed  to  be  essential  to  the  development  of  a successful 
passivation  method. 

7 . 1 Survey  of  Literature  on  Passivation 

There  have  been  a variety  of  technigues  investigated  for  the  growth 
of  native  oxides  or  other  insulating  films  on  gallium  arsenide,  including 
thermal,  anodic,  plasma  oxidation  and  other  methods.  Some  of  these  are  re- 
viewed  by  Schwartz  in  a critical  review  of  GaAs  surface  chemistry  and  by 
Wilmsen,  who  compares  the  growth  of  thermal  oxides  on  the  III-V  compound 
semiconductors  with  anodic  oxides. 

Thermal  oxidation  of  GaAs  has  been  the  subject  of  several  publica- 

tions  , and  one  Ph.D.  research  was  devoted  to  the  thermal  oxidation  of 

3 3 , V*  2 5 

GaAsP.  Murarka  shows  that  amorphous  oxide  films  are  formed  on  GaAs 

heated  in  air  and  that  the  films  have  good  adherence  and  can  be  removed  by 

MCI,  indicating  that  they  might  provide  a useful  mask  for  GaAs  device  tech- 

77 

nology.  Navratil  studied  the  oxidation  of  GaAs  in  dry  oxygen  over  the 


45 


temperature  range  400-530°C  and  also  obtained  amorphous  films.  Later, 

28 

Navratil,  Ohlidal  and  Lukes  proposed  a model  for  the  growth  of  thermal 

oxides  on  GaAs  within  the  temperature  range  480-530°C  and  studied  the  sys- 

28 

tern  using  optical  reflectance  measurements.  In  1962,  Minden  studied  the 
oxidation  of  GaAs  in  pure  oxygen  and  in  water  vapor  containing  oxygen  at 
temperatures  between  600°C  and  900°C,  and  observed  electron  diffraction 

patterns  which  indicated  that  the  films  were  crystalline  and  composed  of 

10 

e-Ga;,03.  Later,  Rubenstein  found  that  the  thermal  oxidation  of  GaAs  in 

pure  oxygen  at  temperatures  between  700°C  and  1130°C  produced  crystalline 

ii 

fl-Ga;,03  and  a small  amount  of  crystalline  GaAsO,,.  Zaininger  and  Revesz 

found  that  GaAs  oxidized  at  660°C  in  dry  oxygen  produced  granular  oxides 

12 

free  of  arsenic.  Also,  Sealy  and  Hemmtnt  concluded  that  some  trace  of 

arsenic  was  present  in  oxide  films  grown  in  oxygen  at  610°C  but  those 

11,34 

grown  at  680°C  contained  no  arsenic.  Finally,  Coerver  found  that  the 
thermal  oxidation  of  GaAsP  in  dry  oxygen  in  the  temperature  range  between 
600°C  and  90CTC  comprised  a mixture  of  crystalline  B-Ga^  and  GaPOi, . Using 
high  energy  resolution  Auger  and  X-ray  photoelectron  spectroscopy  (XPS  or 
ESCA)  techniques,  we  have  shown  (Chapter  V of  this  report)  that  thermal 
oxides  grown  in  oxygen  at  600  °C  were  Ga?03  and  that  they  contained 
elemental  arsenic  throughout  the  bulk  and  a thin  layer  of  As?03  at  the 
oxide-air  interface. 

The  growth  of  anodic  oxides  on  GaAs  has  frequently  been  investigated, 

35-17 

principally  by  Schwartz  at  Bell  Labs  and  by  Hartnagel  and  associates 

at  the  University  of  Newcastle  upon  Tyne. 38-41  Some  of  the  basic  ideas  on 
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the  anodic  oxidation  of  metals  and  semiconductors  were  given  in  1972  by 

4 4 i»  S 

Fehlner  and  Mott.  In  1976,  Weimann  and  Schlapp  reviewed  the  two  tech- 
niques which  give  homogeneous  uniform  films  of  low  conductivity,  constant 
voltage  anodization  in  aqueous  H202  (due  to  Schwartz)  and  anodization  in  a 
mixture  of  glycol  and  water  (due  to  Hartnagel ) . These  electrolytes  were 
also  compared  with  several  non-aqueous  electrolytes  to  determine  the  ex- 
perimental conditions  for  stable  and  uniform  films.  Spitzer,  Schwartz  and 
3e 

Weigle  report  that  anodic  films  grown  in  hot  concentrated  (30%)  H202  can 
be  etched  by  HC1  and  H20  but  that  their  chemical  stability  is  enhanced  by 
heating  in  dry  N2  at  250°C  for  two  hours,  after  which  H20  will  no  longer 
etch  the  oxide.  They  suggest  that  a device  process  for  improved  reliability 
might  utilize  an  oxide  growth-strip-regrowth  cycle  for  removing  surface 
contamination.  Hasegawa,  Forward  and  Hartnagel  investigated  other  aqueous 
electrolytes  than  the  mixture  of  glycol  and  water;  these  included  tartaric 
or  citric  acid  mixed  with  glycol  and  a weak  carboxylic  acid  with  poly- 

38 

hydric  alcohol  in  water.  In  the  second  case,  they  obtained  electrical 
properties  of  the  anodic  films  by  MOS  capacitance  measurements  and  found 
that  the  growth  technique  yields  stable  and  reproducible  layers  with  a 
breakdown  strength  of  5*10f'  volts/cm  and  resistivities  of  1014-1016  ohm-cm. 
Also,  low  temperature  annealing  in  hydrogen  improved  the  interface  proper- 

4? 

ties.  In  1976,  El-Safti,  Weiss  and  Hartnagel  discussed  the  anodization 

43 

of  aluminum  on  GaAs  and  Weiss  and  Hartnagel  showed  that  GaAs  surfaces 
cleaned  in  organic  solvents,  polished  in  bromine-methanol  solutions  and 
etched  in  HC1  have  a residual  thin  layer  of  oxide  on  them.  This  is  con- 
sistent with  the  studies  we  have  performed  using  AFS  to  provide  much  more 
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detail  on  the  composition  of  the  residual  film. 

Schwartz,  Ermanis  and  Brastad  recently  demonstrated  the  anodic 
oxidation  of  GaAs  in  water  of  the  proper  conductivity  and  pH,  and  showed 
that,  for  device  fabrication,  a photoresist  could  be  used  to  define  re- 
stricted areas  for  oxidation  and/or  etching.  Other  electrolytes  studied 
include  a non-aqueous  solution  of  ammonium  pentaborate  by  Arora  and  Bid- 

46  4 7 4 8 

nurkar  and  ammonium  phosphate  in  water  by  Law  and  Lee.  Szpak  re- 
lated the  initial  two-dimensional  phase  of  anodic  oxide  growth  to  a dis- 

4 9 

tribution  of  oxide  nuclei  on  the  surface.  Coleman,  Shaw  and  Dobrott 

obtained  the  oxygen  concentration  profile  in  anodic  films  using  radioactive- 

so 

labeled  oxygen,  while  Guenthner,  Rao,  Thomas  and  Paulson  used  simultaneous 
ion  sputtering  and  Auger  analysis  to  obtain  composition  profiles  for  anodic 
films  grown  on  GaAsP  in  both  aqueous  and  non-aqueous  electrolytes , and 
changes  in  the  profiles  after  annealing. 

There  have  been  few  studies  of  the  growth  of  oxides  in  GaAs  in  an 

SI  6 2 

oxygen  plasma.  Weinreich  used  a technique  similar  to  that  of  Ligenza 
who  studied  the  oxidation  of  silicon  in  a microwave-generated  plasma.  His 

o 

films  were  amorphous  and  reached  thicknesses  of  3500A  in  100  minutes  in  an 

c *1 

oxygen  partial  pressure  of  0.2  Torr.  However,  Sugano  and  Mori  obtained 
polycrystalline  films  grown  on  (100)  and  (111)  GaAs  in  a high  frequency 
oxygen-discharge  plasma,  with  a dependence  on  0^  partial  pressure  and  the 
high  frequency  field.  Composition  profiles  were  made  of  these  films  by 

5 4 5.; 

simultaneous  ion  sputtering  and  Auger  analysis.  Chang  and  Sinha  ob- 
tained amorphous  layers  by  their  plasma  growth  technique  but  they  report 
that  interface  properties  need  to  be  improved  before  they  can  be  used  in 


48 


56 

the  passivation  of  devices.  Chang,  Chang  and  Sheng  recently  used  a new 
approach  for  growing  insulating  layers  on  GaAs;  they  first  deposited  alumi- 
num on  GaAs  then  oxidized  the  aluminum  to  A1203  in  an  oxygen  plasma.  The 
resulting  films  were  analyzed  by  AES  and  transmission  electron  microscopy. 

S7 

Kauffman,  Feldman,  Poate  and  Chang  then  used  helium  ion  backscattering 

and  ion-induced  x-ray  techniques  to  determine  the  composition  profile  but 

found  unresolved  discrepancies  with  the  Auger  data. 

Lukes  studied  the  kinetics  of  the  oxidation  of  gallium  arsenide  in 

so  < 

air  at  room  temperature  and  Loeschke  and  Kuehn  studied  the  oxidation  ki-  1 

netics  and  orientation  dependence  of  the  air  oxidation  of  GaAs  doped  with 

tel lurium  or  tin. 

Two  dielectric  films  which  have  proved  successful  on  silicon,  silicon 
dioxide,  Si02,  and  silicon  nitride,  Si3N4,  have  been  explored  for  the  passi- 

80-71  60 

vation  of  GaAs.  Lorenz  studied  MIS  devices  using  Si 02  as  the  gate 

insulator.  He  observed  considerable  dispersion  in  the  C-V  curves  with 
frequency  and  associated  the  high  surface  state  density  and  distribution  of 

time  constants  to  the  insufficient  purity  of  water  used  in  the  cleaning 

6 1 

operations  prior  to  the  Si02  deposition.  Sato  and  Ikeda  used  both  SiO 

6? 

and  Si02  in  MOS  capacitors.  In  1970,  Kern  and  White  reported  on  the  first 
comprehensive  study  of  silica  films  on  GaAs.  The  Si 02  films  were  deposited 
primarily  by  the  pyrolysis  of  tetraethyl  siloxane  under  various  temperature 
and  pressure  conditions.  C-V  measurements  of  MOS  capacitors  were  made  as  a 
function  of  several  variables  including  the  oxide  deposition  conditions,  a 
post-deposition  heat  treatment,  substrate  dopant,  the  crystallographic  ori- 
entation of  the  GaAs  and  the  chemical  treatment  prior  to  deposition  of  Si 02 . 
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The  lowest  value  of  surface  or  interface  state  density  of  1011  states  cm-2 
eV-1  was  achieved  by  using  (111)  GaAs  plus  a combination  of  specific  chem- 
ical treatments  followed  by  rapid  deposition  of  Si02  at  a low  temperature. 

63 

Weissman  deposited  both  Si02  and  Si  3N()  in  an  R.F.  glow  discharge.  From 
his  C-V  measurements  he  deduced  that  both  types  gave  large  densities  of  in- 
terface states.  The  density  near  the  mid-gap  was  between  4-6'1012  cm-2  eV-1 
and  near  the  band  edges  it  was  estimated  to  be  greater  than  1013  cm-2  eV-1. 

64 

Ito  and  Sakai  fabricated  inversion  channel  GaAs  MISFET's  using 
alloyed  source  and  drain  n+  regions  and  a CVD  double  layer  film  of  A1 ^0 3 
and  SiOv  for  the  gate  oxide,  for  improved  stability.  Semushkina,  Marakhonov 

65 

and  Seisyan  produced  films  of  either  A1203  or  Si02  also  by  CVD  techniques. 
Their  C-V  measurements  indicated  majority  carrier  (n)  depletion  due  to  elec- 
tron capture  by  surface  states  supposedly  introduced  during  the  process  heat 
treatment.  The  passivation  of  microwave  Au-GaAs  Schottky-barrier  diodes  by 

6 fj 

R.F.  sputtered  Si02  films  was  studied  by  Stareev,  Piskorski  and  Majewski . 

They  formed  windows  in  the  oxide  layer  and  vacuum  deposited  gold  to  form  the 

67 

Schottky  barrier  after  cleaning  the  GaAs  surface  by  ion  sputtering.  Messick 
formed  insulating  films  of  silicon  oxynitride  on  LPE  layers  on  GaAs  and  fab- 
ricated MISFET's  which  showed  little  hysteresis  and  mobilities  greater  than 
2000  cm2/volt  sec. 

63 

In  addition  to  the  studies  of  Si3Ni,  films  by  Weissman  and  the  silicon 

67  68  69,70  71 

oxynitride  films  by  Messick,  Sugano  et  al , Ward  and  Seki  et  al. 

investigated  the  passivation  of  gallium  arsenide  with  Si  3N,, . Sugano,  Sakai 

68 

and  Mori  deposited  Si  3N,,  films  on  GaAs  surfaces  at  low  temperatures  using 
an  R.F.  glow  discharge  in  a mixture  of  silane  and  nitrogen  gas.  However, 
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even  though  the  technique  is  promising,  the  films  were  not  dense  enough  to 

6-?  70 

have  sufficient  insulating  properties.  Ward  and  Cooper,  Ward  and  Schwartz 
reported  on  investigations  of  CVD-grown  films  of  Si  3N<4  on  GaAs  including 
their  evaluation  by  admittance-voltage  measurements  of  MIS  capacitors  over 
a wide  frequency  range  (lOKHz  to  50MHz).  They  interpreted  their  results  in 
terms  of  large  densities  of  surface  states  distributed  across  the  interface 
into  the  insulator. 

Insulating  layers  of  A1 203  show  promise  for  good  dielectric  films  on 

68 

GaAs.  Sugano,  Sakai  and  Mori  who  reported  the  growth  of  Si3N4  films  in  an 
R.F.  glow  discharge  also  investigated  ti,a  growth  of  A1203  films  on  GaAs  and 
found  that  their  characteristics  depended  upon  the  growth  rate.  As  mentioned 

64 

earlier,  Ito  and  Sakai  obtained  improved  stability  of  GaAs  MISFET's  by 
using  double  layer  films  of  A1203  and  SiO?.  Semushkina,  Marakhonov  and 

6 5 

Seisyan  observed  that  the  majority  carrier  (n)  concentration  was  depleted 

at  the  surface  for  CVD  A1203  films  as  it  was  for  Si O2  films,  probably  for  the 

7? 

same  reason,  the  substrate  heat  treatment.  Singh  and  Hartnagel  were  able 
to  grow  A1 3 layers  on  GaAs  by  evaporating  aluminum  through  a partial  pres- 
sure of  oxygen  with  a resulting  interface  state  density  that  was  lower  than 

73 

for  native  oxides.  Siekanowicz,  Huan,  Hoffmann,  Jolly  and  Chiang  applied 

aluminum  oxide  passivation  to  electron-beam  semiconductor  diodes  and  Wada, 

74 

Yanagisawa  and  Takanashi  showed  that  heat  treatment  of  A1 -GaAs  Schottky 
barrier  diodes  during  and  after  A1  deposition  improved  their  characteristics 
with  respect  to  surface  states. 
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7.2  Proposed  Method 

Based  upon  the  experience  from  previous  insulating  film  studies, 
as  reported  in  the  literature,  and  upon  the  results  of  our  fundamental  in- 
vestigations of  solid  state  reactions  at  interfaces,  we  propose  the  following 
method  for  the  passivation  of  GaAs  surfaces  for  device  applications. 

The  steps  in  the  process  would  include: 

(1)  The  choice  of  a particular  orientation  of  the  GaAs  substrate  and 
the  use  of  semi -insulating  Cr-doped  substrates  because  of  the 
desirability  of  fabricating  field  effect  devices. 

(2)  The  MBE  growth  of  a semi-insulating  or  undoped  buffer  layer  for 
isolation  from  substrate  defects  and  surface  contamination. 

(3)  The  MBE  growth  of  a doped  layer  which  will  become  the  conducting 
layer  in  field  effect  devices.  The  choice  of  dopant  type  and 
concentration  profile  may  also  be  important  for  device  applica- 
tions. 

(4)  A wet  chemical  cleaning  treatment.  Appendix  A gives  the  proced- 
ure we  have  used  for  cleaning  GaAs  substrates  prior  to  MBE  growth 
or  prior  to  surface  studies  using  AES. 

(5)  Growth  of  a thin  anodic  film  as  an  aqueous  electrolyte  to  displace 
the  interface  and  to  provide  an  initial  film,  e.g.,  Ga?0j. 

(6)  A possible  heat  treatment  in  dry  N?  to  stabilize  the  initial 
film  or  in  H?/N?  mixtures  to  improve  interface  properties. 

(7)  Deposition  of  pure  aluminum  to  reduce  or  partially  reduce  the 
initial  oxide  film  and  to  form  A1 ? O3 - Other  metals  may  also 
prove  to  be  feasible  here  if  the  thermodynamics  is  favorable. 
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(3)  Growth  of  thicker  insulating  layers  of  Al;,0-j  by  a suitable 
deposition  process,  for  example,  reactive  sputtering. 

(9)  The  possible  low  temperature  growth  of  Si02  or  other  insulator 
on  the  A1?03  film. 

(10)  A post-deposition  heat  treatment. 
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VIII.  SURFACE/ I NTERFAC E STATES  ON  ELEMENTAL  AND  COMPOUND  SEMICONDUCTORS : 

(A.  Craig  Kenton) 

8.1  I_n  t roduction 

Thus  far  we  have  considered  problems  associated  with  the  physics 
and  chemistry  of  surfaces  or  interfaces  on  semiconductors  from  a macrosco- 
pic point  of  view.  In  our  studies  during  the  last  year  we  have  focused  on 
the  microscopic  or  atomistic  approach,  especially  in  discussing  the  nature 
of  interface  reactions  which  will  depend  upon  the  local  atomic  arrangements 
and  chemical  bonds.  The  influence  of  these  local  effects  on  the  electronic 
properties  of  the  semiconductor  is  reflected  in  the  energies  and  densities 
of  surface  or  interface  states. 

A study  was  initiated  under  the  program  of  theoretical  techniques  for 
describing  the  localized  electronic  structure  at  the  atomic  sites  of  a quasi- 
molecular  complex  comprising  substrate  and  interface  atoms.  Three  particu- 
lar theoretical  techniques  appear  promising.  These  are:  (1)  molecular  clus- 
ter calculations;  (2)  the  pseudopotential  method;  and  (3)  the  generalized 
Wannier  function  method.  Section  8.5  of  this  chapter  reviews  these  local- 
ized theoretical  approaches. 

Of  the  several  surface  analytical  techniques  which  probe  the  electronic 
structure  of  solids,  photoemission  and  energy  loss  spectroscopy  have  pro- 
vided the  most  lucid  observations  of  intrinsic  surface  or  interface  states. 
Thus,  data  from  photoemission  and  energy  loss  spectroscopy  studies  will  be 
useful  for  comparison  with  the  results  of  the  theoretical  studies.  These 
experimental  techniques  and  the  interpretations  of  surface/interface  states 
for  silicon  and  gallium  arsenide  are  reviewed  in  Sections  8.2  and  8.3. 

The  techniques  of  LEED  and  RHEED  provide  information  on  the  crystal- 
lographic order  and  lattice  damage  of  a clean  surface  or  surface/adsorbate 
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system.  The  theoretical  analysis  of  sir.h  data,  however,  is  formidable  due 
to  complications  arising  from  surface  relaxation  and/or  reconstruction,  mul- 
tiple scattering  processes  and  many  body  effects.  Much  progress,  however, 
has  been  made  recently  in  this  area  yielding  useful  model  geometries  that 
can  be  inserted  into  the  electronic  structure  calculations.  Models  for  the 
atomic  geometry  of  gallium  arsenide  surfaces,  based  primarily  on  LEED  data, 
are  discussed  in  Section  8.4. 

8 . 2 Surface/Interface  States  and  Photoemissiori  Spectroscopy 

In  principle,  the  photoemission  experiment  involves  using  electro- 
magnetic radiation  in  the  ultraviolet  to  X-ray  range  as  an  exciting  probe 
to  eject  electrons  out  of  a solid.  The  ejection  of  electrons  will  only 
occur  if  the  exciting  photon  energy  is  sufficiently  high.  Once  outside  the 
solid,  electrons  may  be  analyzed  in  terms  of  their  kinetic  energy  and  emis- 
sion angle.  Photoemission  probes  the  intrinsic  electronic  properties  of 
the  sample  yielding  a spectrum  which  will,  in  general,  reflect  both  bulk 
and  surface  contributions.  Surface  sensitivity  of  photoemission  spectra  on 
surfaces  can  arise  from  two  theoretical  conditions  resulting  in  a localized 
electronic  wavefunction:  (1)  trapping  at  a surface  state  localized  in  a 
band  gap;  and  (2)  strongly  damped  final  state  wavefunctions ,75 ,7S  The  last 
condition  arises  from  the  inelastic  electron-electron  interactions  which  re- 

O 

suit  in  electron  scattering  lengths  of  ' 5 A,  yielding  emitted  photoelec- 
trons in  the  energy  range  20  - 100  eV.  This  condition  is  responsible  for 
the  surface  sensitivity  of  any  electron  spectroscopy,  e.g.,  Auger  Electron 
Spectroscopy  (AES)  or  Low  Energy  Electron  Loss  Spectroscopy  (LELS),  in  this 
range . 

8.2.1  Studies  on  Silicon  (Si)  Surface/Interface  States 

Bardeen77  first  postulated  the  existence  of  surface 
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electronic  states  on  semiconductor  surfaces  to  explain  band  bending  and 
Fermi  level  pinning  at  a metal -semiconductor  contact  (Schottky  barrier). 
Allen  and  Gobeli™  observed  that  the  photoelectric  threshold  and  work  func- 
tion of  Si (111)  remained  constant  as  the  Fermi  level  was  swept  through  the 
bulk  band  gap  by  varying  the  impurity  concentration,  i.e.,  the  Fermi  level 
remained  pinned  at  the  surface  due  to  the  existence  of  surface  states. 

They  deduced  that  the  surface  states  existed  in  two  bands,  one  above  and 
one  below  the  Fermi  level;  however,  since  the  photon  energy  was  only  about 
1 eV  above  the  valence  band,  the  full  energy  distribution  of  the  surface 
states  was  not  obtained.  Chiarotti,  et  al  .7’  observed  optical  absorption 
for  photon  energies  less  than  the  bulk  band  gap  by  maximizing  the  surface 
absorption  sensitivity.  After  the  cleaved  surface  was  studied,  it  was  then 
exposed  to  oxygen  and  restudied.  The  differences  in  the  clean  and  oxidized 
absorption,  shown  in  Figure  17,  was  assigned  to  surface  states.  Direct  pho- 
toemission evidence  of  filled  surface  states  on  Si  (111)  was  first  obtained 
by  Eastman  and  Grobman80  and  by  Wagner  and  Spicer81  for  photon  energies 
more  than  5 eV  above  the  threshold.  The  results  clearly  showed  a shoulder 


Valence  Band  Maximum 

Figure  17.  Model  of  Si  Surface  States  as  Determined  by 

Optical  Absorption  (after  Chiarotti  et  al.  )• 


near  the  top  of  the  valence  band  at  the  surface  which  disappeared  after  low 
exposure  to  air  or  residual  gas.  Small  differences  in  the  location  of  peaks 
have  been  attributed  to  the  well  known"''  (2x1)  reconstructed  surface  of 
Si(lll)  after  cleavage  and  the  (7  x 7)  annealed  structure.85  Figure  18 
illustrates  the  experimental  spectra  of  several  authors.  Spicer87  suggests 
one  method  of  identifying  filled  surface  states  is  by  varying  the  photon 
energy  and  observing  the  shifting  of  bulk  peaks.  This  is  presently  most 
conveniently  done  by  using  continuously  tunable  synchrotron  radiation.  The 
establishment  of  the  pinning  of  the  Fermi  level70  and  verification  by  more 
thorough  photoemission  studies80  illustrates  how  band  bending  affects  sili- 
con surface  states. 


bp  low  ZjUvj  r /p  fcV) 

Figure  18.  Comparison  of  Photoemission  Energy  E1str1 bution  Spectra 
of  a)  Freshly  Cleaved  S f (111)  2x1  and  b)  Annealed 
(>  850°C)  Si(lll)  7x7  Surfaces  by  Several  Authors00'®* 
(after  Feuerbacher  and  Willis70). 

Figure  19  shows  the  shift  of  bulk  states  with  respect  to  the  surface  Fermi 
level  and  the  relatively  minor  shift  of  the  surface  states.  An  experimental 
approach  of  separating  bulk  and  surface  photoemission  data  by  averaging  over 
the  different  surface  structures  and  subtracting  the  average  from  the  indi- 
vidual spectra  yielding  the  surface  contribution  is  illustrated  in  Figure 

8 9-91 

20.  Justification  of  this  procedure  is  that  the  averaged  spectrum 
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Figure  18.  (a)  Energy  Distribution  Curves  for  Silicon.  The  high- 

energy  peak  is  due  to  electrons  in  surface  states.  The 
pinning  effect  causes  a 0.2  eV  movement  of  this  peak 
and  a small  shift  in  the  work  function  edge  (left). 

(b)  Band  Diagram  for  n++  (left)  and  p++  (right)  Silicon 
Illustrating  the  Pinning  Effect  (after  Spicer”^). 


Figure  20.  (a)  Photoemission  energy 

distribution  from  Si(lll)  7x7  with 
the  estimated  contribution  from  secondary 
electrons  removed,  (b)  Comparison  of 
"averaged"  photoemission  curve  with  the 
density  of  states  calculated  by  Kane92. 

(c)  Surface  transition  density  of  states 
obtained  by  Rowe  and  I bac hB 3 by  subtracting 
the  assumed  bulk  signal  (curve  b)  from 
the  measured  energy  distribution  (curve  a). 


CNERCY.(.V) 
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is  similar  to  a theoretical  calculation  of  the  density  of  states.9*  This 
appears  most  useful  in  identifying  surface  resonances,  i.e.  surface  states 
which  are  degenerate  with  hulk  states. 

Theoretical  progress  on  silicon  has  been  achieved  by  several  authors  ,9  ,_9 '' 
notably  by  Applebaum  and  Hamann, 9e"9B  where  they  have  correlated  one 

of  the  surface  states  in  the  band  gap  with  "dangling  bonds"  of  the  surface 
atoms.  Taking  into  account  self-consistency  and  surface  relaxation,  two 
surface  states  are  found. The  (2x1)  reconstruction  of  freshly 
cleaved  Si(lll)  has  been  theoretically  studied  by  several  authors  .9  9 ’ 1 0 0-1 ''  1 

The  use  of  angle-resolved  photoemission,  i.e.  the  energy  analysis  of 
ejected  electrons  in  terms  of  their  emission  angle  suggests  new  interpre- 
tations of  the  bonding  mechanisms  of  surface  species.  Although  many  experi- 
mental and  theoretical  difficulties  still  exist,  results  for  Si (111)  indi- 
cate appreciable  mixing  of  the  dangling  bonds  and  back  bonds  . 1 0 4 0 9 

In  addition  to  conventional  energy  analysis  of  photoemission  data,  the 
technique  of  photoemission  yield  spectroscopy  has  been  used  to  probe  un- 
filled bulk  and  surface  states.  For  Si(lll),  the  work  of  Koma  and  Ludeke1 0 1 
and  Bauer,  et  a 1 . 1 1 ® i s cited. 

A widely  used  method  for  the  characterization  of  surface  states  on 
semiconductors  is  by  chemisorption  to  form  surface  bonds.  The  difference 
in  the  clean  surface  spectra  and  the  chemisorbed  spectra  are  then  compared. 

In  the  case  of  silicon,  it  is  found  that  the  dangling  bonds  react  with 
several  active  gases,  preferentially  oxygen,  resulting  in  the  disappearance 
of  the  observed  intrinsic  surface  states.  On  n-type  samples,  oxidation 
removes  the  silicon  surface  states  and  destroys  the  Fermi  level  pinning. 

The  pinning  is  not  totally  destroyed  on  p-type  silicon  even  though  the 
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surface  state  peak  vanishes  from  the  photoemission  spectra.  The  remaining 
small  density  of  surface  states  now  become  interface  states  resulting  in  the 
partial  pinning  on  p-type  samples.88*111  This  dramatic  effect  (see  Figure 
21) is  the  primary  reason  that  the  Si/SiO^  interface  is  a passive  structure — 
providing  reliable  exterior  surfaces  for  the  countless  devices  utilizing 
silicon  technology.  Interface  states  have  been  studied  theoretically  and 
experimentally  on  the  formation  of  the  monohydride  phase  [Si(lll)  7x7+ 

Si  (1 1 1 ) : H]1 1 1 and  the  trihydride  phase  [Si(lll)  lxl-+  Si  ( 1 1 I ) : Si  J ‘ ' 
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Figure  21.  The  Effect  of  Oxygen  on  the  Si  Surface  States.  The  peak 
at  about  -1.0  eV  (0  L)  is  due  to  surface  states,  which 
have  almost  entirely  disappeared  after  101  or  10I(L 
exposure.  (1  Langmuir,  L = 10"6torr-sec .)  (after  Spicer11  *) . 
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of  silicon.  The  correlation  between  atomic  geometry  and  electronic  struc- 
ture for  this  sytem  has  proved  fruitful.  Very  recently,  the  atomic  position 
of  Cl  chemisorbed  on  Si(lll)  has  been  studied  using  angle-resolved  photo- 
emission with  synchrotron  radiation.115  The  electronic  structure  has  been 
correlated  with  the  atomic  structure  by  the  pseudopotential  method. 

The  study  of  silicon  surface/interface  states  provides  the  concepts 
by  which  we  will  extend  our  knowledge  to  more  complicated  systems,  for  ex- 
ample, the  compound  semiconductors.  The  various  techniques  of  photoemission 
continue  to  play  a crucial  role. 

3.2.2  Studies  on  Gallium  Arsenide  (GaAs)Surface/Interface  States 

Studies  concerning  surface  and  interface  states  on  GaAs  have 
unquestionably  practical  significance  --  the  surface  passivation  of  gallium 
arsenide,  for  example,  is  one  area  receiving  much  current  attention.  Due  to 
the  fact  that  GaAs  is  a compound  semiconductor,  one  would  expect  it  to  be 
much  more  complicated  than  the  elemental  semiconductors.  Early  work  function 
measurements  indicated  this  fact  by  the  absence  of  Fermi  level  pinning  by 
surface  states. 115A  large  band  bending  for  n-type  samples  and  the  absence 
of  band  bending  on  p-type  samples  further  documented  the  complexity.117 
Conclusions  reached  were  that  no  arsenic-derived  surface  states  existed  with- 
in the  bulk  band  gap  and  that  an  unfilled  band  of  gallium-derived  surface 
states  existed  in  the  gap.  Eastman  and  Grobman,1 11  using  photoemission,  con- 
tradicted these  results  by  observing  a filled  band  of  surface  states  for 
GaAs(llO)  lying  in  the  band  gap.  It  was  later  shown  that  steps  on  badly 
cleaved  surfaces  could  induce  extrinsic  surface  states  in  the  gap.11 9,1 20  Sub- 
sequent photoemission  results  to  date  have  generally  agreed  on  the  absence 
of  any  filled  intrinsic  surface  states  within  the  band  gap. 121-1  The 


61 


existence  of  the  unfilled  states  in  the  band  qap,  however,  has  been  the 

1 24 

subject  of  much  controversy.  Huijser  and  Van  Laar  concluded,  by  contact 
potential  difference  measurements,  that  no  empty  surface  state  band  exists 
in  the  band  gap  of  the  GaAs(llO)  cleaved  surface.  This  result  contradicted 

120-2  125-130  1211 31-4 

several  experimental  ’ and  theoretical  ’ studies  but  now 

1 20  1 3 5-7 

has  been  confirmed  by  photoemission  experiments.  ’ The  present 

status  of  the  GaAs(llO)  surface  is  perhaps  best  illustrated  by  the  modified 
GSCH  (Gregory-Spicer-Ciraci-Harrison)  model  shown  in  Figure  22. 

In  agreement  with  the  absence  of  filled  surface  states  in  the  bind  gap, 

l ?0  12?  130 

no  Fermi  level  pinning  is  observed  for  p-type  GaAs(lin)  ’ ’ supporting 

li  c 

earlier  studies.  Until  recently,  photoemission  and  other  techniques  fa- 
vored the  pinning  of  the  Fermi  level  for  n-type  GaAs  at  the  midgap. 


Figure  22.  Schematic  GSCH  Model  of  the  GaAs(llO)  Surface  Showing 
the  Energy  Position  of  Filled  and  Empty  Surface  States 
and  Their  Association  with  the  Ga  and  As  Surface  Atoms 
(after  Lindau  et  al.121  ). 
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Huijser  and  Van  Laai^2"  persisted  with  the  earlier  views  and  maintained  that 
their  GaAs(llO)  cleaves  on  n-type  material  showed  no  pinning  of  the  Fermi 
level.  Furthermore,  Chye,  et  a 1 ,m  have  shown  the  absence  of  pinninq  on 
both  n-  and  p-type  gallium  antimonide  (GaSb).  Pinning  was  shown  to  occur 
if  a small  amount  of  oxygen  was  admitted  to  the  system  or  if  a Schottky 
barrier  was  formed.  Gudat,  et  al  ,137  have  also  found  III -V  n-type  compounds 
with  an  absence  of  pinning.  These  facts  motivated  Spicer  and  workers1*’1''1'1 
to  study  the  question  of  n-type  GaAs(llO)  in  detail.  In  addition  to  surface 
defects,  strain,  or  steps  induced  by  cleaving,  they  found  that  a small  amount 
of  contamination  on  a good  cleave  would  stabilize  the  Fermi  level  around  the 
midgap.  The  substantiated  result  is  that  Fermi  level  pinning  does  not  occur 
on  cleaved  n-type  GaAs(llO)  due  to  intrinsic  surface  states.  Pinning,  how- 
ever, does  occur  on  these  surfaces  due  to  extrinsic  causes.  GaAs(llO),  there- 
fore, does  not  possess  intrinsic  surface  states  of  any  type  in  the  bulk  band 
gap.  This  is  in  agreement  with  Figure  22. 

It  is  noteworthy  to  state  that  the  above  revelations  came  about  essen- 
tially by  detailed  studies  on  interface  states,  primarily  concerning  the 
uptake  of  oxygen  on  GaAs(llO),  using  synchrotron  radiation  as  the  exciting 
source  for  photoemission,136’1'*'''2  Pianetta,  et  a 1 .l4?  suggest  that  the  slight 
contamination  which  caused  erroneous  results  on  good  cleaves  was  caused  by 
the  excitation  of  oxygen  due  to  the  ionization  gauge  and  ion  pump  operating 
in  the  ultra-high  vacuum  (UHV)  environment.  Previous  studies  showed  an  oxy- 
gen uptake  rate  of  10^  - 10^  L (1L  = 1 Langmuir  = 10-fl  Torr-sec)  for  a half- 
monolayer and  a subsequent  slow  sorption  rate.  Taking  the  precautions  of 
valving  off  the  ion  pump  and  using  a cold-cathode  (Redhead)  gauge,  Pianetta, 
et  al.147  have  shown  a marked  exposure  increase  (It/  - 10^  L)  required  for 
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half-monolayer  coverage.  Also,  the  monitoring  of  the  As-3d  and  Ga-3d  photo- 

emission  core  levels  has  shown  that  the  As  peak  shifts  but  that  the  Ga  peak 

1 2 

remains  fixed  for  exposures  to  1 0 L.  After  this  exposure,  the  Ga-3d  peak 
only  broadens  slightly,  as  shown  in  Figure  23.  Previous  experiments  have 
shown  that  the  As  peak  shifts  and  that  the  Ga  peak  remains  fixed  up  to 
in6  L exposure  but  that  the  Ga  peak  shifts  at  higher  exposures,  as  shown 
in  Figure  24. 


OXIDATION  or  n-TYf*t  GoAt(llO)  AT  100  «V 
Go  3d 


Figure  23.  Photoemission  Spectra  of  Clean  and  Oxygen  Exposed 

GaAs(llO)  Surfaces  at  a Photon  Energy  of  100  eV.  Note 
the  2.9  eV  shifted  As  3d  peak  and  the  high  oxygen  exposure 
101 2L,  without  any  shift  of  the  Ga-3d  peak  (after 
Pianetta  et  al J4* ) . 
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Figure  24.  Photoeml ssion  Spectra  of  GaAs(llO)  Surfaces  Exposed  to 
Activated  Oxygen.  Note  the  appearance  of  a second 
shifted  (AE  = 4.5  eV)  As-3d  peak  and  a 1.0  eV  shifted 
Ga  3d  peak  (after  Lindau  et  al.123). 


The  later  data.  Figure  23,  suggest  that  a major  charge  transfer  takes  place 

from  arsenic  to  oxygen,  illustrating  that  oxygen  is  preferentially  bonded 

to  the  As  surface  atoms  and  not  to  the  Ga  atoms  on  the  truly  clean  surface. 

1 2 

The  fact  that  the  Ga  peak  remains  fixed  for  the  long  exposure  of  10  L 
(equivalent  to  atmospheric  pressure  of  0^  for  20  minutes)  suggests  that 
the  bonds  between  the  surface  atoms  and  the  bulk  are  not  broken.  So,  for 
a short  time,  a freshly  cleaved  GaAs  specimen  at  room  temperature  and  atmos- 
pheric pressure  should  possess  stable  back  bonds. 123,135  The  earlier  results, 
Figure  24,  showing  the  shifting  of  the  Ga-3d  peak,  are  attributed  to  oxygen 
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broken  up  by  a hot  filament  ionization  gauge,  attacking  the  bonds  binding 

1'tS 

the  surface  atoms  to  the  rest  of  the  crystal . The  movement  of  the  Fermi 
level  in  these  studies  has  been  monitored  during  the  oxygen  exposure  and 
is  shown  in  Figure  2b  for  a total  of  10  cleaves.  The  figure  shows  that  sur- 
face perfection  and  minute  oxygen  contamination  dramatically  affects  the 
initial  stages  of  oxidation.  The  inset  in  the  figure  shows  that  the  valence 
band  photoemission  data  from  samples  beginning  with  the  Fermi  level  unpinned 
and  moving  to  the  midgap  have  much  sharper  features  than  the  initially  pinned- 
samples1 23  ’ 1 1,f)  Lindau,  et  al.l41and  Pianetta,  et  al.140  have  recently  studied 
in  detail  the  valence  band  structure,  using  photoemission,  of  cleaved  and 
oxygen  exposed  GaAs(llO).  In  addition  to  the  results  that  have  been  men- 
tioned, they  emphasize  that  the  chemisorption  of  small  amounts  of  oxygen 
(<  10  L)  produces  long  range  ordering  effects  which  sharpen  the  intrinsic 
valence  photoemission  spectra.  This  is  attributed  to  rearrangement  within 
the  unit  cell  resulting  in  the  relieving  of  surface  strain. 


Figure  25.  Fermi  Level  Pinning  on  GaAs(llO)  for  Various  Samples 

and  Cleavage  Qualities  as  a Function  of  Oxygen  Exposure. 
Samples  1 and  II  are  supplied  by  Laser  Diode  Corp. 

(doping  3 x 1017  cm-3)  andsample  III  by  MCP  Electronics, 
Ltd.  (doping  1.7  x 10*B  cm”3).  The  insert  shows 
representative  photoemission  spectra  for  samples  with 
pinned  and  unpinned  Fermi  level  (after  Lindau  et  al.123). 
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In  the  study  of  the  metal -semiconductor  interface,  the  Bardeen1'*1'  model , 
which  correlates  Fermi  level  pinning  with  the  density  of  intrinsic  surface 


states,  has  been  used  extensively.147’1 14,1  However,  the  fact  that  a metal  is 
added  to  a surface  has  raised  many  theoretical  questions  concerning  the 


continued  existence  of  surface  states.  1 49-1  54  In  the  case  of  cleaved  GaAs(llO), 


the  previous  discussion  has  established  that  for  both  n-  and  p-type  material, 
the  surface  Fermi  level  is  unpinned.  After  a metal  is  added,  however,  the 
Fermi  level  is  found  to  be  pinned  at  the  interface.1 16  This  suggests  that 
the  metal  induces  new  interface  states,  now  commonly  referred  to  as  metal 
induced  gap  states  (MIGS).  Heine149  suggests  that  the  new  states  are  due 
to  the  tails  of  the  metallic  wavefunction  decaying  into  the  semiconductor. 
Inkson191  suggests  that  the  states  are  due  to  the  narrowing  of  the  semicon- 
ductor gap  at  the  interface.  Phil  1 ips1  ,;i<  proposed  that  the  states  ar  e due 
to  elementary  excitations  and  chemical  bonding  that  pin  the  Fermi  level. 
Perhaps  the  latest  theory  has  been  presented  by  Mele  and  Joannopoul os.1 9 9 
It  correlates  the  Fermi  level  pinning  to  the  formation  of  a dipole  layer 
which  is  due  to  the  filling  of  dangling  bond  states.  Thus,  the  intrinsic 
surface  states  are  still  shown  to  play  a vital  role  in  the  formation  of  the 
Schottky  barrier.  They  point  out  that  their  theory  is  presently  valid  only 
for  tetrahedral ly  coordinated  semiconductors  and  suggest  methods  of  exten- 


sion. 


The  GaAs/Cs  interface  has  been  studied  by  photoemission  by  Spicer  and 
co-workers1 30 »’ 95-7  in  several  instances.  However,  the  experimental  studies 
that  were  undertaken  are  essentially  invalid  for  a true  metal -semiconductor 
junction  due  to  their  later  conclusions  regarding  post-cleavage  contamina- 
tion of  their  surfaces.  A detailed  study  of  the  variation  in  barrier  height 
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of  Schottky  barriers  as  a function  of  the  interface  oxide  thickness  appears 
to  be  very  feasible  using  photoemission  techniques. 


•4 

In  summary,  intrinsic  surface  states  on  cleaved  GaAs(llO)  now  are  in 
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Oxidation  studies  on  GaAs(llO)  have  been  characterized  in  terms  of  the  ini- 
tial and  the  long  term  uptake  and  are  correlated  well  with  the  pinning  of 
the  Fermi  level  and  the  location  of  the  bonding  sites.  Tunable  synchrotron 
radiation  provides  a clean  and  wide  range,  stable  energy  source  for  photo- 
emission  experiments  under  excellent  vacuum  conditions. 

In  the  fabrication  and  passivation  of  elemental  silicon  for  electronic 
devices,  Fermi  level  pinning  is  substantially  reduced  when  the  surface  un- 
dergoes oxidation.  As  a result,  the  devices  are  useful  and  reliable.  The 
opposite  behavior  is  encountered  for  the  oxidation  of  GaAs(llO),  resulting 
in  the  Fermi  level  being  pinned  in  the  band  gap.  It  is,  therefore,  not  sur- 
prising that  a good  passivation  scheme  has  not  yet  been  developed. 

8.3  S ur fa ce/Interface  States  a nd  El ectron  Energy  Loss  Spectroscopy 
Interest  in  Electron  Energy  Loss  Spectroscopy  (ELS)  as  a surface 
analytical  technique  has  come  into  use  within  the  last  decade  due  to  the 
ability  to  prepare  reproducible  surfaces  in  an  ultra-high  vacuum  environment. 
Experiments  using  ELS  have  concerned  collective  vibrational  states  (phonons), 
normal  vibrational  states,  and  electronic  states  of  clean  and  gas-covered 
surfaces.  This  discussion  is  related  to  the  study  of  the  excitation  of 
electronic  surface  transitions  on  semiconductor  surfaces. 

lSft 

Most  theoretical  results  of  ELS  are  based  on  the  dielectric  theory.' 

A high  density  of  surface  states  may  cause  the  dielectric  constant  for  a 
surface  to  differ  from  that  of  the  bulk.  The  surface  is  then  modelled  as 
a thin  layer  possessing  a surface  dielectric  constant  which  depends  on  the 
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distance  normal  to  the  vacuum-solid  interface.  Anisotropic  effects,  how- 
ever, have  been  experimentally  observed.1™ 

Typically,  the  energy-loss  spectra  are  measured  with  a conventional 
Auger  single-pass  cylindrical  mirror  analyzer  (CMA)  having  an  integral 
coaxial  electron  gun.  The  primary  electron  beam  strikes  the  sample  at  nor- 
mal incidence  with  energies  in  the  50-1000  eV  range.  At  these  low  energies 
(typically  < 500  eV)  the  escape  depth  of  the  secondary  electrons,  due  to 
inelastic  scattering  processes,  is  small,  resulting  in  surface  sensitive 
spectra.  When  ELS  is  performed  in  this  way,  it  is  referred  to  as  Low- 
Energy  Electron  Loss  Spectroscopy  (LLLS)  and  one  operates  in  the  reflection 
mode  of  ELS.  The  energy  resolution  required  for  the  CMA  is  100  meV  based 
on  the  magnitude  of  expected  energy  losses.  To  obtain  this  resolution, 
modulation  voltages  are  kept  small  (0.4  - 1.0  volts,  peak-to-peak) . To 

increase  the  sensitivity,  the  second  harmonic  of  the  modulation  voltage  is 

d^N 

detected  to  obtain  second  derivative,  — r-  , spectra,  using  standard  lock-in 

dE^ 

techniques.  It  should  be  noted  that  the  differentiation  annihilates  un- 
structured loss  intensity  and  overestimates  sharp  structures  of  low  intensity. 

Eroitzheim,  et  al  .1W  studied  cleaved  and  annealed  surfaces  of  Si(lll) 
for  the  first  time,  calculated  optical  absorption  data,  and  compared  it  to 
existing  data  .la»16?  They  noticed  the  vanishing  of  intrinsic  surface  tran- 
sitions after  annealing.  Ibach,  et  al.161  studied  the  LELS  spectra  of  Si(lll) 
and  Si (100)  2 x 1 as  a function  of  oxygen  coverage  up  to  a monolayer.  They 
noticed  the  disappearance  of  surface  transitions  of  Si  for  the  clean  surfaces, 
the  appearance  of  surface  transitions  related  to  oxygen,  the  decrease  in  in- 
tensity of  bulk  and  the  surface  plasmons,  and  a splitting  of  the  surface 
plasmon,  all  as  a function  of  increasing  oxygen  coverage. 
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Further  LELS  studies  on  Si1  f'R'1  70 and  Ge1  f;n r4_1  ■>  ‘ /!  show  comparable  re- 


sults. Of  particular  mention  are  the  combined  photoemission  and  LELS  experi- 
ments by  Rowe,  et  al.170  performed  on  Si (111)  enabling  an  energy  level  model 


Froitzheim  and  Ibach177  have  studied  LELS  for  the  cleaved  and  oxygen 


covered  GaAs(llO)  surface  and  have  not  found  any  structure  that  can  be  asso- 


ciated with  surface  transitions.  This  has  been  attributed  to  low  step  den- 
sities and  the  lack  of  significant  wavefunction  overlap  between  As  and  Ga 

surface  states.  These  results  are  in  agreement  with  the  latest  photoemission 

173  at 

experiments  performed  by  Lindau,  et  al.  Stanford. 

In  direct  contradiction  to  these  experiments,  Ludeke,  et  al . 1 c4-r-  1U 

have  observed  transitions  involving  the  narrow  filled  Ga-3d  core  levels  as 

initial  states  and  empty  surface  levels  or  conduction  band  levels  as  final 

states  on  GaAs(lOO),  (111),  and  (110)  surfaces.  In  Figure  26,  their  data 

show  a peak  at  21.3  eV  for  clean  GaAs(lOO)  and  Ga-stabil ized  (111)  surfaces 

which  is  associated  with  excitations  from  the  Ga-3d^,j  core  level  into  empty 

surface  states.  A weaker  peak  at  20  eV  is  a similar  transition  from  the 

Ga-3dr/0  core  level.  This  structure,  attributed  to  spin-orbit  splitting, 
b/ z 

disappears  as  the  surface  becomes  As-stabilized  and  a strong  peak  appears 
at  1.7  eV  associated  with  transitions  from  filled  As  surface  states  to  empty 
conduction  band  states. 


Room  temperature  oxygen  adsorption  on  various  surfaces  was  determined 

from  the  ratio  of  the  peak-to-peak  amplitudes  of  the  0(510  eV)  and  Ga(1070 

eV)  signals  in  Auger  dN/dE  spectra  as  a function  of  exposure- ’ 74 ’ 1 7 Mono- 

layer  coverage,  as  indicated  by  saturation  in  the  0/Ga  Auger  ratio,  is  com- 

g 

plete  for  exposures  of  the  order  of  10  L.  The  Ga-3d  to  surface  state  peaks 
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Figure  26.  Second-Derivative  Energy- 
Loss  Spectra  of  GaAs:  (a)  (100) 
Surface,  Ga  Stabilized;  b)  (TIT) 
Surface,  Ga  Stabilized;  c)  (111) 
Surface,  As  Stabilized;  d)  (III) 
Surface,  Sputtered  (after  Ludeke 
and  Esaki15  4 ) . 
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are  observed  to  diminish  in  intensity  corresponding  to  the  oxygen  adsorp- 
tion. This  effect  is  interpreted  as  saturation  of  the  unfilled  Ga  dangling 
bonds  of  the  clean  surface.  Simultaneously  with  the  disappearance  of  the 
above  peaks,  new  peaks  appear  at  20.9  and  23.2  eV,  which  are  believed  to 
be  due  to  excitations  from  the  Ga-3d  core  levels  into  empty  oxide-related 
interface  states.  No  changes  are  observed  in  the  loss  peaks  due  to  the  ex- 
citation of  As-3d  core  levels.  It  is  thus  concl uded1 68 74-5 that  in  the 
oxidation  of  (100)  and  (111)  surfaces,  as  well  as  other  GaAs  surfaces  in 
general,  the  oxygen  is  chemically  bonded  to  Ga  atoms  and  that  the  As  atoms 
not  directly  involved  in  the  oxidation  process. 

'ne  different  GaAs  surfaces  were  prepared  either  by  Ar+  sputtering  and 
I-  ; r by  molecular  beam  epitaxial  techniques.  The  As  and  Ga  surface 


71 


compositions  of  (100)  surfaces  were  readily  changed  by  treatments  involving 
cooling  a cleaned  (100)  surface  below  200  °C  in  an  arsenic  vapor  of  10  ^ 
Torr  (As  rich  surface),  heating  in  ultra-high  vacuum  to  400  °C  (As  stabi- 
lized), or  heating  t.o  about  500  °0  (As  depleted).  Reflection  high  energy 
electron  diffraction  (RHEED)  patterns  establish  the  surface  crystallography 
of  these  three  different  (100)  surfaces17'  Low  energy  electron  diffraction 
(LEED)  could  also  be  used  to  determine  the  surface  reconstruction  in  each 
case.  When  the  surfaces  are  prepared  by  Ar+  sputtering  and  annealing,  one 
would  have  to  monitor  the  ratio  of  the  Ga-to-As  low  energy  (MW)  Auger  peaks 
as  a measure  of  the  surface  composition  and  the  ratio  of  the  Ga-to-As  high 
energy  (LMM)  Auger  peaks  as  a measure  of  the  bulk  composition.175"177 

The  previous  conclusion  concerning  the  site  of  oxygen  chemisorption 
has  been  the  subject  of  much  controversy178’ 179  since  photoemission  results 
clearly  indicate  that  the  As  atoms  are  the  initial  bonding  site173  Lapeyre 
and  Anderson180  have  presented  strong  evidence  that  excitonic  effects  are 
involved  in  the  core  1 evel -surface  state  transition.  They  have  proposed  a 
direct  recombination  process,  in  which  the  core  hole  and  electron  recombine 
and  transfer  their  total  energy  to  valence  electrons.  This  would  require  a 
strong  localization  of  the  initial  electron-hole  pair  suggesting  the  exci- 
tonic nature  of  the  empty  surface  state.  An  alternative  model  involving 
internal  photon  creation  by  the  relaxation  of  the  core  level  has  also  been 
proposed  by  J.  C.  Phillips.  These  observations  suggest  that  the  energy 
locations  of  the  empty  surface  states  of  GaAs(llO)  must  be  mod  if  ied.181 
When  this  is  done,  the  results  of  LELS  are  in  some  accord  with  the  recent 
photoemission  experiments.173  Photoemission  results  indicate  GaAs(llO) 
possesses  no  intrinsic  surface  states  (filled  or  unfilled)  in  the  band  gap 
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and  oxygen  initially  bonds  to  the  As  surface  atoms.  Mele  and  Joannopoulos:  3 
present  a theoretical  argument  using  an  empirical  tight-binding  model  which 
indicates  that  oxygen  bonds  mostly  to  the  As  surface  atoms  with  some  degree 
of  mixing  to  the  Ga  species.  Detailed  points  of  the  apparent  controversy 
between  the  photoemission  and  LELS  results  are  still  unresolved.178’179 

L ELS  studies  on  AlAs  have  also  been  performed.164  This  is  an  important 

O 

compound  due  to  the  fact  that  the  bulk  lattice  spacing  (5.6622  A)  is  very 

o 

close  to  that  of  GaAs  (5.6532  A).KE  AlAs  overlayers  can  also  easily  be 
grown  on  GaAs  substrates  by  MBE,  resulting  in  minimal  interface  distortion 
due  to  lattice  mismatching.177 

With  respect  to  metal  overlayers  on  GaAs,  Rowe186  has  used  a combination 
of  ultraviolet  photoemission  spectroscopy  (UPS)  and  electron  energy-loss 
spectroscopy  to  study  the  effect  of  metal  layers  on  (111)  GaAs  surfaces. 

UPS  of  direct  bulk  transitions  gives  information  on  changes  in  the  surface 
potential  on  adsorption  and  energy-loss  measurements  show  the  resultant 
changes  in  the  occupied  and  unoccupied  intrinsic  surface  states.  In  addi- 
tion, Gudat,  Eastman,  and  Freeouf167  report  on  empty  surface  states  on  the 
(110)  surfaces  of  several  III-V  semiconductors,  including  GaAs,  and  the 
effect  of  metal  overlayers  on  these  states,  using  photoemission.  In  both 
of  these  studies  there  is  a correlation  between  the  surface  state  energies 
and  the  position  of  the  Fermi  level  (or  value  of  the  surface  potential)  at 
the  metal -semiconductor  interface.  Thus,  it  is  concluded,  as  a number  of 
model  studies  suggest,  (see,  for  example,  the  recent  paper  by  Borrego, 
Gutmann,  and  Ashok)188  that  the  intrinsic  surface  or  interface  states  still 
play  an  important  role  in  determining  Schottky  barrier  heights. 

The  relatively  new  method  of  LELS  yields  another  surface  analytical 
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technique  which  is  important  in  the  analysis  of  the  microstructure  of  solid 
surfaces.  The  contradictions  that  have  arisen,  in  comparison  to  other  experi- 
ments, have  required  a closer  scrutiny  of  the  effects  of  some  standard  experi- 
mental procedures  and  have  motivated  better  explanations  of  the  spectra  from 
theorists.  LELS,  still  in  its  infancy,  will  most  likely  provide  useful  ex- 
perimental inputs  into  the  complex  theories  of  low  energy  electron  scattering. 

8. 4 Surface  Atomic  Geometry  of  GaAs  (110) 

In  every  useful  theoretical  model  which  seeks  to  calculate  local 
properties  of  matter,  a detailed  knowledge  of  the  geometrical  configuration 
of  atomic  species  at  the  surface  or  interface  region  is  required.  The 
analysis  of  elastic  low-energy  electron  diffraction  (ELEED)  data  has  pro- 
vided most  of  this  knowledge  to  date. 

Surface  atomic  geometry,  in  principle,  can  be  determined  from  the  analy- 
sis of  LEED  intensities.189  Surface  sensitivity  is  achieved  because  of  the 
strong  interaction  of  low  energy  electrons  with  the  first  few  atomic  layers 
of  a solid.  This  introduces  the  problem  of  multiple  scattering  effects  which 
makes  dynamical  LEED  calculations  complex  and  lengthy.  In  practice,  there 
exists  a multiplicity  of  possible  geometries  which  will  yield  a particular 
LEED  structure.  The  usual  procedure  is  then  to  calculate  the  LEED  inten- 
sities from  probable  geometries  and  compare  the  results  to  experimental  in- 
tensities. The  atomic  geometry  which  best  fits  the  experimental  data  is 
then  selected.  Even  though  the  GaAs  (110)  surface  has  the  ideal  1X1  LEED 
structure,  two  factors  suggest  that  the  atomic  arrangement  of  the  surface 
is  different  than  that  of  the  bulk.190  The  (hk)  and  (hk)  beams  exhibit  ex- 
treme asymmetry  in  intensity  and  the  (10)  and  (10)  beams  possess  strong 
intensity.  If  the  surface  was  an  ideal  termination  of  the  bulk  zincblende 
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lattice,  the  (10)  and  (10)  beams  would  be  very  weak. 

The  following  discussion  will  outline  the  model  for  the  GaAs  (110)  re- 
constructed surface  proposed  by  Duke,  et  al  .l91  Assuming  that  covalent  bonds 
are  relatively  incompressible,  one  can  model  the  surfaces  of  covalent  semi- 
conductors by  altering  bond  angles  while  keeping  bond  lengths  fixed.  In 
support  of  this  assumption,  the  reduction  in  surface  Madelung  energy  and 
the  breaking  of  surface  covalent  bonds  suggest  a reduction  of  charge  transfer 
in  the  surface  region  compared  to  the  bulk.  For  GaAs  (110)  1X1,  these  assump- 
tions lead  to  the  hypothesis192*193  that  the  As  atoms  rotate  outward  from  the 
surface  tov/ard  a prismatic  coordination  and  that  the  Ga  atoms  rotate  inwards 
toward  a three-fold  planar  (sp)  coordination.  The  rearrangement  occurs  only 
in  the  first  layer,  as  shown  in  Figure  27,  and  is  completely  characterized  by 
a single  parameter,  9y,  the  tilt  angle.  This  is  the  angle  between  the  (110) 
surface  and  the  plane,  S,  passing  through  the  nearest  neighbor  surface  Ga  and 
As  atoms.  The  figure  shows  the  ideal  1X1  structure  (9-j.  = 0)  alongside  the 
reconstructed  1X1  model.  Some  interatomic  distances  are  shown  along  with 
the  1X1  surface  net  in  the  top  view.  In  the  side  view,  the  decrease  in  the 
average  interlayer  spacing  between  the  surface  and  second  layer  is  denoted 
as  L (which  is  directly  related  to  9j  in  this  model).  This  rearrangement 
preserves  not  only  the  bond  lengths  but  also  the  periodicity  parallel  to  the 
surface.  This  simplification  is  very  important  as  several  theoretical  me- 
thods (e.g.,  the  pseudopotential  method)  rely  on  surface  translation  symme- 
try in  order  to  perform  "layer"  calculations.  Indeed,  Duke,  et  al  .191  found 
that  their  exact  multiple-scattering  computer  programs194  required  too  much 
memory  for  zincblende  or  wurtzite  cleavage  face  calculations.  New  programs 
were  developed,  based  on  the  layer  iteration  method,  where  individual  layer 
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scatterin']  amplitudes,  reflection  matrices,  and  transmission  matrices  are 
first  calculated.19*  1,15 


ATOMIC  STRUCTURE  OF  GoAs(IIQ)  Ixl  SURFACF 
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Figure  27.  Atomic  Structure  of  Ideal  and  Relaxed  GaAs(llO)  1 x 1 

Surfaces.  The  ideal  surface  is  obtained  by  truncating 
the  bulk  lattice.  In  the  relaxation  model,  Ga  atoms 
move  into  and  As  atoms  move  out  of  the  crystal  (after 
Duke  et  al.  )• 


From  their 
with  experiment 
with  the  larger 


calculations,  Duke,  et  al  .lql  determined  the  best  agreement 
occurred  for  the  range  of  the  tilt  angle  27°  < oT  < 34.3°, 
angles  being  statistically  preferred. 
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Duke,  et  al  .ni  also  suggest  a correlation  between  the  local  electronic 
structure  and  the  atomic  geometry  of  the  surface  region  mainly  due  to  the 
character  of  the  chemical  bonding,  i.e.,  the  percent  ionicity  of  covalency 
of  the  bond.  Pandey,  et  al  ,n7  have  used  the  previously  discussed  geometric 
model  in  tight-binding  calculations  to  generate  the  local  density  of  states 
for  the  GaAs  (110)  1X1  surface.  Studying  two  values  of  the  tilt  angle, 
f>T  = 19°  and  o^.  = 35",  they  concluded  that  0^  = 19°  compared  more  favorably 
with  various  photoemission  experiments.  It  is  of  interest  to  note  that 
their  local  density  of  states  curves,  shown  in  Figure  28,  show  a sharp  peak 
just  below  the  conduction  band  minimum  (E  ) for  their  best  choice  of  0^  = 19". 


Figure  28.  A Comparison  of  the  Bulk 
Density  of  States  with  the  Local  Densities 
of  States  at  the  Surface  for  the  Ideal, 
the  P.elaxed(l  °° ) , and  the  Fully  Relaxed 
(35°)  Models  of  GaAs(llO)  lxl  Surface 
(after  Pandey  et  al.1'*7). 
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This  is  in  disagreement  with  the  latest  photoemission  data,19”  which  shows  no 
intrinsic  surface  states  in  the  bulk  band  gap.  However,  the  fully  relaxed 


surface,  o 
s tci 1 0 , in 


= 35°,  shows  the  movement  upward  and  broadening  of  this  surface 
rssmsn  t with  p ho  to  9m  i s s i o n ddtd  3P.d  t ho  LEED  calculations  o f 


Duke,  et  al . 

Very  recently,  Tong,  et  a l.199  have  compared  their  experimental  LEED 
intensity  data  with  curves  calculated  by  the  combined  space  method  ( CSM ) ? 0 0 
This  dynamical  method  for  calculating  LEED  intensities  is  specifically  de- 
signed for  closely  spaced  noncoplanar  surface  layers.  Extending  the  pre- 
viously described  model,  they  have  found  that  intensity  curves  are  sensitive 
not  only  to  the  tilt  angle  but  also  to  the  vertical  height  of  the  surface 

i 

atoms  from  the  second  layer.  This  corresponds  to  viewing  A as  an  indepen- 
dent parameter  in  the  previous  model.  The  best  agreement  with  experiment, 
determined  from  11  relaxed  structures,  was  found  to  have  a 9y  = 27°  and  a 

o o 

surface  layer  separation  distance  of  1.452  A (a  = 0.233  A).  This  distance 
corresponds  to  a slight  compression  of  the  surface  layer  Ga  and  As  back  bonds. 

The  bond  lengths  of  the  surface  atoms  to  one  another  remain  fixed  for  this 
model . 

For  comparison  to  previous  calculations,  Tong,  et  al.  studied  the  geo- 
metries with  tilt  angles  6y  = 19°  and  0y  = 34.8°.  They  found  considerable 
disagreement  with  the  published  results,  attributing  it  to  an  inaccurate 
method  in  calculating  the  multiple  scattering  effects. 

The  basic  model  that  was  described  by  Duke,  et  al . and  extended  by 
Tong,  et  al . remains  intact  to  date  and  provides  the  best  information  known 
on  the  structure  of  the  GaAs  (110)  1X1  surface.  With  all  the  bond  lengths 
fixed,  it  has  been  used  for  various  theoretical  calculations;  in  particular, 
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the  pseudopotential  simulation  of  the  GaAs  (110)  1X1  surface.701 

The  importance  of  the  geometrical  structure  cannot  be  underestimated 
in  theoretical  calculations.  Better  theoretical  models  and  subsequent  cal- 
culations of  LEED  intensity  curves  appear  to  be  the  most  useful  endeavor* 
especially  to  understand  the  differences  in  the  results  of  the  present  cal- 
culations. The  strong  correlation  between  surface  atomic  geometry  and  local 
electronic  structure  are  increasingly  becoming  a method  by  which  assumed 
model  geometries  are  conceived  and  tested  by  LEED  calculations.  Recently, 
Lindau,  et  al  .2B2  have  suggested  that  the  surface  electronic  structure  may 
provide  more  detailed  information  on  rearrangement  of  surfaces,  especially 
GaAs  (110),  than  LEED.  This  assertion  needs  to  be  substantiated  by  quanti- 
tative calculations. 

8.5  Theoretical  Modeling  of  Localized  Electronic  States 
3.5.1  Hartree-Fock  Molecular  Cluster  Calculations 

The  evolution  of  experimental  and  theoretical  surface  analy- 
tical techniques  over  the  past  twenty  years  has  resulted  in  significant  in- 
creases in  resolution,  sensitivity,  and  accuracy.  In  conjunction  with  this 
trend,  the  concepts  incorporated  into  the  band  theory  of  solids,  primarily 
developed  to  understand  gross  and  bulk  features  of  crystalline  matter,  have 
necessarily  been  modified  in  order  to  accomodate  microscopic  phenomena  such 
as  impurities,  lattice  vacancies,  etc.  Furthermore,  macroscopic  phenomena 
such  as  crystal  surfaces,  have  required  extensive  reformulation  of  several 
aspects  of  band  theory.  The  band  theory  approach  has  been  and  will  continue 
to  be  highly  successful.  On  the  other  hand,  approaches  dealing  with  the  as- 
semblage of  atoms  to  form  a solid  from  "first  principles"  are  conceptually 
and  intuitively  attractive.  The  ultimate  success  of  both  techniques  would 
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be  their  ability  to  simultaneously  model  a particular  experimental  system. 

The  advent  of  high  speed  computers  has  moved  theoretically-established  first 
principle  calculations  of  practical  systems  into  the  realm  of  possibility. 

The  most  important  and  useful  model  that  has  been  developed  for  elec- 
tronic structure  calculations  is  the  Hartree-Fock  (Hf)  method , which  is 
intimately  connected  with  the  one-electron  theory  of  matter.  P.  W.  Anderson"'01 
has  stated  that  the  main  usefulness  of  the  fundamental  many-body  theory  is 
predicting  the  range  of  application  and  flexibility  of  the  Hartree-Fock 
theory.  In  calculating  the  total  energy  of  an  atom  in  its  ground  state, 
the  error  in  the  actual  experimental  energy  and  the  calculated  Hartree-Fock 
energy  is  defined  as  the  correlation  energy.  The  percentage  error  in  the 
correlation  energy  steadily  decreases  from  '1%  for  He(Z=2)  to  0.1%  for 
Ar(Z=18)  and  appears  to  continue  its  downward  trend  for  heavier  elements. 20,4 
The  fact  that  such  a small  correction  is  all  that  is  necessary  to  yield  ex- 
perimental results  shows  that  the  Hartree-Fock  method  is  a very  good  approxi- 
mation to  the  many-electron  problem. 

In  extending  the  Hartree-Fock  method  to  molecules,  the  major  pitfall 
is  the  destruction  of  the  spherical  symmetry  of  the  system.  Solutions  of 
the  central -field  problem  (i.e.  atomic  orbitals)  then  become  invalid  and 
molecular  orbitals  must  be  introduced.  The  general  procedure  for  handling 
polyatomic  problems  can  be  traced  back  to  Roothaan's  method,  also  referred 
to  as  the  analytic  Hartree-Fock  method. 204  In  this  method,  electron  wave- 
functions  are  modeled  as  a linear  combination  of  a finite  number  of  suitably 
chosen  bas^i_s  functions . The  problem  then  becomes  one  of  algebraically  de- 
termining the  coefficients  of  the  basis  functions  in  each  wavefunction  ex- 
pansion. The  procedure  commonly  used  for  solving  the  Hartree-Fock  equations 
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is  the  method  of  the  self-consistent  field  (SCF).  This  iterative  procedure 
can  be  carried  out  to  any  predetermined  degree  of  accuracy  by  comparing 
successive  calculations  of  the  set  of  coefficients  for  the  basis  functions. 

The  choice  of  a suitable  set  of  basis  functions  to  model  the  one-electron 
wavefunction  has  been  the  subject  of  extensive  investigation. 

Perhaps  the  most  well  known  basis  function  occurs  in  the  "minimal -basis 
set,"  and  is  referred  to  as  a Slater-type  orbital  (STO): 

n-1  -ar  „ > 

r e 

where  Y?m(e,<(>)  is  a spherical  harmonic.  For  a central-field  problem,  such 
as  a single  atom  or  a molecule  having  infinite  separation  between  atoms, 
modeling  of  a single  one-electron  wavefunction  by  a STO  is  very  appropriate.704 
In  order  to  model  a real  molecule,  which  possesses  finite  internuclear 
distances,  the  parameter  a (the  orbital  exponent)  has  been  introduced.  A 
single  one-electron  wavefunction,  in  this  case,  is  still  modeled  by  one  STO 
but  now  a is  optimized  by  minimizing  the  energy  of  the  system. 

A logical  extension  of  this  method  is  the  "double-zeta"  basis  set. 

Here,  a single  wavefunction  is  modeled  as  a linear  combination  of  two  STO's. 

If  reasonable,  but  different,  orbital  exponents  are  selected  for  the  two 
STO's,  a simpler  linear  ootimization  of  the  coefficients  of  the  STO's  is 
all  that  is  necessary  to  give  comparable  results  to  a minimal  basis  set  cal- 
culation. 

The  basis  functions  that  have  been  discussed  can  be  extended  to  include 
"polarization  functions"  (basis  functions  with  higher  angular  quantum  num- 
bers), which  allow  charge  buildup  in  the  bonding  region.200  It  is  now  clear 
that  the  flexibility  introduced  by  this  refinement  is  important  in  geometri- 
cal structure  predictions. 
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Furthermore,  completely  general  functions  can  be  used  for  a single 
orbital  allowing  greater  flexibility,  but  with  a considerable  increase  in 
computational  difficulties.207 

STO's  have  been  used  with  great  success  in  atomic  and  molecular  cal- 
culations. In  Lrie  molecular  problem,  however,  one  is  faced  with  the  addi- 
tional complication  of  assigning  the  one-electron  orbitals  to  particular 
geometric  centers  associated  with  the  equilibrium  sites  of  the  atoms.  The 
analytic  Hartree-Fock  procedure  then  requires  one  to  calculate  two-electron 
integrals  over  all  orbitals  in  the  system.208  This  presents  a tremendous 
calculational , but  tractable,  problem  for  functions  like  the  STO's  in  cal- 
culating the  many-center  two-electron  integrals  for  a molecular  problem. 

To  facilitate  this  difficult  problem  of  integral  computation,  S.  F.  Boys209 

p m n *ar2 

suggested  the  use  of  the  cartesian  Gaussian  functions  of  the  form:  X‘Y  Z e 
The  distinct  advantage  of  the  Gaussian  functions  as  a basis  set  relies  on 
the  elegant  theorem  that  a product  of  two  Gaussian  functions  is  also  a 
Gaussian.  As  a result,  the  most  complicated  four-center  two-electron  inte- 
grals can  be  reduced  immediately  to  two-center  cases  which  can  rapidly  be 
evaluated  by  simple  formulas.  However,  Gaussian  functions  are  less  approp- 
riate than  functions  such  as  STO's  for  representing  one-electron  orbitals. 

As  a result,  a linear  combination  of  Gaussian-type  orbitals  (GTO's)  must  be 
used  for  obtaining  comparable  accuracy.  Typically,  3-8  GTO's  are  used  for 
a particular  one-electron  wavefunction  expansion.  The  use  of  GTO's,  then, 
creates  a substantial  increase  in  the  number  of  integrals  to  be  evaluated. 

For  two-electron  integrals,  the  increase  varies  as  the  fourth  power  of  the 
number  of  GTO's  in  an  expansion.  Furthermore,  the  final  calculation  of  the 
electronic  wavefunction  is  proportional  to  the  number  of  integrals.  Even 
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i ith  the  increase  in  the  number  of  calculations,  it  generally  pays  to  use 
GTO's  rather  than  STO's  for  efficient  molecular  calculations  due  to  the 
significant  decrease  in  the  time  per  calculation. 

Satisfied  with  using  a linear  combination  of  GTO's  as  a basis  set  for 
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cients  to  describe  a one-electron  orbital.  One  method  of  doing  this  is  by 
expanding  a minimum  basis  STO  into  a linear  combination  of  K GTO's.710 
The  resulting  basis  sets  are  denoted  by  STO-KG  (Slater-type  orbital  at  the 
K-Gaussian  level).  These  are  documented  for  K=2  through  6 for  several  ele- 
ments.711 As  K increases,  the  STO-KG  will  approach  the  true  STO.  ST0-3G 
and  ST0-4G  yield  comparable  results  to  a full  STO  calculation. 

To  combat  the  large  number  of  Gaussian  basis  functions  that  accumulate 
as  one  studies  a heavier  atom  or  a larger  molecular  cluster,  several  "con- 
traction" schemes  have  been  introduced.71 7 Contracted  functions  are  formed 
by  taking  fixed  linear  combinations  of  the  Gaussian  basis  functions,  now 
called  primitives.  Care  must  be  exercised  in  constructing  compact  and  flex- 
ible contractions.  Generally,  inner  shell  orbitals  are  contracted  to  some 
degree  while  leaving  outermost  orbitals  uncontracted.  This  procedure  is 
called  a segmented  contraction  scheme.  A general  contraction  scheme  has 
been  used  by  Raffenetti713  which  reduces  all  core  orbitals  to  a single  func- 
tion. The  valence  orbitals  are  modeled  by  a linear  combination  of  core  pri- 
mitives and  the  innermost  valence  function.  The  remainder  of  the  valence 
orbitals  are  left  uncontracted. 

Optimized  atomic  Gaussian  basis  set  expansions  are  presently  available 
through  xenon  (Z=54).714  Segmented  contraction  schemes  which  give  good  re- 
sults have  been  established  through  copper  (Z=29).717 
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^ One-electron  wavefunctions , in  general,  possess  a spatial  and  spin  co- 

l ordinate  dependence.  If  there  are  no  constraints  placed  on  the  total  Hartree- 

Fock  wavefunction  (other  than  orthonormality)  and  the  electron  wavefunction 
is  denoted  by  the  antisymmetric  combination  of  a spatial  wavefunction  times 
a spin  eigenfunction,  the  version  of  the  theory  is  referred  to  as  unrestricted 
Hartree-Fock  (UHF).  If  the  wavefunctions  are  constrained  to  be  the  same  for 
different  spins,  i.e.  spin  is  neglected,  the  version  of  the  theory  is  called 
restricted  Hartree-Fock  (RHF).  In  cases  where  a system  has  unpaired  elec- 
trons, the  UHF  theory  usually  leads  to  an  energy  lowering.  For  larger  atoms 
or  molecules,  the  UHF  method  has  been  shown  to  properly  model  the  contraction 
and  localization  of  higher  angular  momentum  electrons.  In  particular,  for 
| transition  metal  hydrides,  UHF  calculations21 5 show  the  proper  localized 

behavior  of  d-electrons,  while  other  methods  such  as  RHF216  and  SCF-Xa-SW217 
(Self-Consistent  Field  X-alpha  scattered  wave)  calculations  show  considerable 
I;  s-p-d  mixing. 

Efficient  computer  programs,  developed  by  Dr.  Isiah  Shavitt  and  his 
group  at  Battelle  Columbus  Laboratories,  are  presently  available  at  Georgia 
Tech  for  calculating  the  complete  electronic  structure  of  a molecular  cluster 
of  atoms.  The  code  is  designed  to  perform  RHF,  UHF,  and  MCSCF218,219  (mul- 
tiple configuration  SCF)  calculations  using  Gaussian  functions  as  a basis 
set.  Some  of  the  features  of  the  program  are:  (1)  the  ability  to  internally 
generate  Gaussian  expansions  of  STO's;  (2)  the  ability  to  flexibly  input  any 
desired  segmented  contraction  scheme  of  given  primitive  Gaussians;  (3)  the 
exact  treatment  of  the  exchange  interaction;  (4)  the  performance  of  charge- 
conserving  integral  approximations;220  (5)  the  efficient  calculation  of  two- 
electron  integrals;220  and  (6)  a contour  and  perspective  granhical  package 
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for  visualizing  the  results.  In  conjunction  with  the  above,  there  is  an 
option  also  available  which  allows  replacement  of  the  core  orbitals  with 
an  ab  initio  effective  core  potential.  This  technique,  developed  and  im- 
plemented primarily  by  L.  Kahn,  enables  one  to  consider  only  detailed  cal- 
culations involving  the  valence  electrons,  significantly  reducing  computa- 
tion time.221 *22M  Using  this  method,  feasible  self-consistent  calculations 
of  valence  orbital  wavefunctions  and  energies  can  be  made  for  heavier  atoms 
and  larger  molecular  clusters. 

The  basis  of  the  molecular  cluster  approach  in  modeling  surface/inter- 
face electronic  states  at  a semiconductor  surface  relies  on  the  assumption 
that  the  electronic  structure  can  be  reasonably  described  by  the  local 
atomic  arrangement.  As  the  cluster  size  is  increased,  the  results  should 
approach  that  of  a semi-infinite  solid  as  described  by  band  structure  cal- 
culations. Unfortunately,  practical  limitations  on  computation  time  restrict 
the  size  of  a cluster  greatly.  Finite  size  clusters  add  the  additional  com- 
plication that  important  differences  may  be  observed  between  the  electronic 
structure  of  a cluster  and  a semi-infinite  solid.  This  is  primarily  due 
to  the  fact  that  most  atoms  are  at,  or  are  very  close  to,  the  surface  of  a 
feasible  cluster.  One  particular  method  that  has  been  used  for  silicon  to 
reduce  effects  of  the  cluster  edges  is  to  introduce  hydrogen  atoms  to  tie 

up  the  "dangling  bonds" . 22''  «22r'  The  hydrogen  atom  is  expected  to  saturate 

3 

the  bonds  for  silicon  since  its  electronegativity  matches  that  of  the  sp 
hybridized  bonding  orbitals  and,  thus,  preventing  any  unwanted  charge  trans- 
fer at  the  outermost  cluster  sites.  Another  method  of  reducing  unwanted 
cluster  edge  effects  has  recently  been  used  by  Nishida.227  He  points  out 
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that  using  H atoms  for  saturating  semiconductor  bonds  is  not  general  and 

is  not  satisfactory  for  germanium  and  other  semiconductor  surfaces.  In  his 

extended  Hiickel  theory  (EHT)229’229  calculations,  he  has  chosen  an  appro- 

3 

priate  basis  set  and  Hamiltonian  to  exclude  the  surface  sp  orbitals  of 
silicon  and  germanium  clusters.  Hydrogen  chemisorption  on  Si(lll)  and 
Ge(lll)  surfaces  as  well  as  the  clean  GaAs(llO)  ideal  and  relaxed  surface 
was  studied.  The  results  obtained  by  EHT  show  evidence  of  surface  and  inter 
face  electronic  states  due  to  charge  localization.  The  major  drawback  of 
EHT  calculations  is  that  it  is  a non-SCF  procedure.  Nevertheless,  it  has 
proven  the  feasibility  of  performing  cluster  calculations  on  semiconductor 
surfaces. 

There  are  several  methods  by  which  the  results  of  cluster  studies  may 
be  applied  to  experimental  and  other  theoretical  studies. 

The  comparison  of  energy  levels  of  the  molecular  cluster  to  band  struc- 
ture results  has  been  the  predominantly  used  procedure.  One  can  determine 
the  band  gap  and  valence-band  width  as  a function  of  cluster  size.  Figure 
29  shows  this  comparison  for  the  EHT  cluster  calculations  of  silicon.  Fur- 
thermore, the  grouping  of  energy  levels  suggest  the  coagulation  of  bands. 
Figure  29  shows  this  grouping  for  silicon.  For  adsorption  studies  on  metals 
such  as  CO  on  Ni,  Figure  30  shows  how  the  density  of  energy  levels  is  quali- 
tatively compared  to  photoemission  band  structure  results.  Comparison  with 
gas  phase  energy  spectra,  of  CO  in  this  case,  shows  the  non-trivial  identi- 
fication and  resequencing  of  energy  levels  upon  adsorption.230 

By  analyzing  the  molecular  wavefunctions  resulting  from  a cluster  cal- 
culation, one  may  calculate  the  local  density  of  states  (LDOS).  This  may 
better  be  compared  to  surface  sensitive  photoemission  results  and  other 
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Figure  29.  Energy  Level  Diagrams  for  Ideal  Si(lll)  Clusters 
as  a Function  of  Cluster  Size:  (a)  SI  1 g ; (b)  Si 2 g 
(c)  S i 4 s ; and  (d)  Si 5 s - The  result  for  the  in- 
finite cluster  is  shown  in  (e)  and  originates  from 
a band  structure  calculation.  Levels  labelled  by 
D correspond  to  sp3  dangling  bond  orbitals  (after 
Nishida227). 


Figure  30.  Comparison  of  Theoretical  Energy  Levels  for  the 
C0(N1)s  Cluster  and  the  Experimental  Photo- 
emission  Spectrum  of  CO  Chemisorbed  on  Ni(100). 
The  inset  shows  the  model  cluster  (after  Batra 
and  Bagus230). 
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theoretical  models.  Figure  31  shows  the  EHT  LDOS  calculations  for  a GaAs(llO) 
cluster  which  is  in  agreement  with  pseudopotential  layer  calculations231 
and  tight-binding  calculations,232  but  not  in  present  agreement  with  pre- 
viously discussed  recent  photoemission  results  in  Section  8.2.2.  Further 
analysis  of  the  molecular  wavefunctions  will  give  information  on  their  angu- 
lar dependence  and  allow  comparison  to  angle-resolved  photoemission  studies. 


Figure  31.  Extended  Hlickel  Theory  Calculations  for  Ideal  (solid  lines) 
and  Relaxed  (dashed  lines)  GaAs(llO)  Surfaces  Modeled  by  a 
GamAsju  Cluster.  First  layer  (a),  second  layer  (b),  and 
third  layer  (c),  LDOS's  are  shown.  In  (d),  the  total  density 
of  states  is  shown,  with  prominent  surface  states  indicated 
by  arrows  (after  Nishida227). 
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A very  important  application  of  the  cluster  method  is  predicting  struc- 
tural geometries  of  a system.210  Applications  of  this  type  that  have  re- 
cently begun  to  appear  in  the  literature  are  primarily  concerned  with  the 
chemisorption  of  a single  adatom  onto  a surface.  A study  of  H chemisorp- 
tion on  a NiO  surface233  has  been  carried  out  using  RHF,  UHF,  and  the 
generalized  valence  band  (GVB)234  methods.  Figure  32  shows  the  total  energy 
of  the  chemisorption  system  versus  the  distance  the  H atom  is  from  the  sur- 
face plane  for  different  adsorption  sites.  The  minima  of  the  potential 
curves  yield  the  binding  energy  and  equilibrium  distance  for  the  H-0  or 
H-Ni  bond.  A similar  study  of  the  chemisorption  of  0 on  the  Li (100)  surface 
has  been  performed  using  RHF.235  Three  sites  were  considered  for  the 


Figure  32.  Potential  Energy  Curves  for  a H Atom  Approaching  a Surface 
of  NiO.  The  upper  curves  correspond  to  H chemisorption 
on  the  0 sites  of  NiO;  the  lower  curves  correspond  to  H 
chemisorption  on  the  N1  sites  of  NiO  (after  Bobrowltz 
and  Goddard234). 
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sorption,  as  shown  in  Figure  33A.  Binding  energy  versus  surface  distance 
results  were  obtained  for  different  cluster  sizes,  as  shown  in  Figure  33B, 

C and  D.  The  bridge  site  for  this  system  was  determined  to  yield  the  more 
stable  configuration,  thus  predicting  a suitable  geometry  for  subsequent  study. 
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Figure  33.  Chemisorption  of  0 on  the  Li(100)  Surface.  A shows  the 
three  high-symmetry  surface  adsorption  sites;  B,  C,  and 
D show  the  oxygen  binding  energy  as  a function  of  distance 
from  the  surface  plane  for  the  on-top,  central,  and  bridge 
site,  respectively  (after  Hermann  and  Bagus235). 
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A thorough  study  of  the  HF  method,  especially  concerning  selection  of  basis 
functions  and  cluster  convergence  properties,  has  recently  been  done  for  H 
chemisorption  on  the  Be(0001 ) surface.236  It  was  found  that  even  a 22  atom 
Be  cluster  did  not  have  properties  of  a bulk  metal.  They  established  that 
the  interaction  of  a H adatom  with  the  substrate  is  predominantly  determined 
by  nearest  Be  neighbors.  Properties  associated  with  chemisorption  are, 
therefore,  expected  to  converge  rapidly  with  respect  to  cluster  size.  Pre- 
dictions on  the  stable  geometrical  location  of  the  adatom  were  carried  out, 
as  previously  described,  by  studying  the  binding  energies. 

As  a result  of  the  feasibility  of  cluster  calculations,  RHF  and  UHF 
studies  were  initiated  under  our  ONR-sponsored  program  on  model  semiconduc- 
tor surfaces  using  the  previously  mentioned  computer  codes.  The  following 
oath  of  approach  has  been  established.  The  compound  semiconductor  GaAs  has 
been  selected  primarily  because  of  current  fundamental  interest  and  the 
existing  difficulties  in  passivating  the  surface  for  device  applications. 
Initial  calculations  are  presently  being  carried  out  on  various  atoms  of  in- 
terest, e.g.  Ga,  As,  0,  A1 , Ta,  and  H,  to  determine  suitable  basis  sets  and 
ah  initio  core  potentials.  Full  optimized  GTO  basis  sets  have  been  obtained 
for  elements  up  to  Xe(Z=54).237  The  next  calculations  will  involve  clusters 
of  Ga  and  As  atoms  arranged  in  various  crystallographic  orientations,  i.e. 
GaAs (110),  (100),  (111)  and  ( T T T ) along  with  attempts  to  reduce  cluster  edge 
effects.  Vacancies  in  the  lattice,  for  example,  resulting  from  arsenic  de- 
pletion during  a heat  treatment,  can  be  considered  by  removing  strategic 
atoms  from  the  cluster.  Subsequent  calculations  would  identify  the  effect 


of  chemisorbed  atoms,  e.g.  Ga  0 , A1  As  , A1  0 , Ta  0 , H,  etc.,  represent- 
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ing  surface  "quasi-molecular"  clusters  from  the  passivation  film  or  interface. 
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In  accomplishing  this,  the  particular  angular  orientation  of  the  electronic 
structure  of  the  added  atoms  must  be  correlated  with  the  underlying  crystal 
structure.  These  will  be  assumed,  mostly  based  on  previous  experimental 
(AES,  LEED,  XPS,  UPS,  etc.)  results,  some  of  which  have  been  described  in 
this  report.  Comparison  of  these  calculations  should  show  detailed  inter- 
face electronic  structure  differences  due  to  changes  in  bonding  configura- 
tions at  the  surface  of  the  initial  cluster.  The  results  will  be  compared 
to  existing  atomic  energy  level  data  and,  in  particular,  photoemission  data. 
Based  on  these  comparisons,  coupled  with  binding  energy  calculations  for 
chemisorbed  clusters,  predictions  can  be  made  on  the  goemetrical  arrangements 
of  the  surface  and  interfacial  layers.  In  line  with  the  prediction  by  Lindau 
(see  Section  8.4)  that  surface  electronic  structure  may  provide  better  in- 
formation on  the  geometrical  rearrangement  of  surfaces,  the  HF  cluster  mo- 
del may  reveal  an  important  theoretical  link  by  which  detailed  rearrangements 
may  be  assumed  and  efficiently  compared  with  experimental  surface/interface 
data. 

The  HE  molecular  cluster  model  provides  a "first  principles"  evaluation 
of  electronic  structure.  The  use  of  this  method  for  the  heavy  atoms,  such 
as  Ga  and  As,  will  require  developing  optimum  small  clusters  for  modeling 
crystal  surfaces  in  order  to  achieve  a practical  computation  time.  The  mo- 
del is  expected  to  yield  relevant  information  on  the  local  electronic  en- 
vironment of  the  interfacial  species  and,  hence,  enable  choices  of  an  opti- 
mum geometric  and  electronic  configuration.  For  passivation  purposes,  the 
optimum  configuration  would  seek  to  reduce  interface  electronic  states. 
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8.5.2  Localized  Electronic  Structure  Calculations  using  Wannier 
Functions 

In  studyino  the  energy  spectrum  in  an  excited  el0  tron  in 
ideal  crystal,  Gregory  Wannier  constructed  "orthogonal  atomic  functions," 
now  called  Wannier  functions.238  The  motivation  for  this  was  that  atomic 
functions,  rather  than  delocalized  Bloch  functions,  had  been  shown  to  be  more 
effective  for  studying  the  behavior  of  electrons  in  crystals.238  The  problem 
with  atomic  functions,  however,  is  that  they  are  not  orthogonal.  We  postu- 
late a function,  an(r),  such  that  a Bloch  function  in  the  nth  band  is  given 
by: 


1 r 


yielding  art  inversion  formula 


It  can  easily  be  shown,  using  the  orthogonality  of  the  Bloch  functions,  that 
the  Wannier  functions,  a n ( v. ) , are  orthogonal  for  different  bands  and  for  dif- 
ferent sites.  Phenomena  involving  localized  electronic  levels  are  thus 
conveniently  represented  in  terms  of  Wannier  functions.  Surface  states  have 
been  studied  utilizing  these  methods.240  Theoretical  progress  has  been  made, 
notably  by  Kohn241  and  Cloizeau.242  The  properties  of  Wannier  functions  in 
relation  to  energy  bands  are  firmly  established  in  these  studies.  Useful 
quantitative  calculations  had  not  appeared  until  1972  when  Kohn243  adapted 
an  ah  initio  variational  technique  — calculating  Wannier  functions  for 
several  bands  of  the  diamond  structure.  Once  Wannier  functions  are  found. 
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Bloch  functions,  energy  eigenvalues,  and  moments  of  the  density  of  states 
are  easily  calculated.  This  method,  then,  provides  an  alternative  approach 
to  conventional  bulk  band  structure  calculations. 

A milestone  occurred  when  Kohn  and  Onffroy  (1973)?l*‘*  introduced  the 
concept  of  Generalized  Wannier  Functions  (GWF)  in  studying  a one-dimensional 
periodic  lattice  with  a point  defect.  They  have  shown  that  the  GWF's  have 
the  same  localized  properties  as  the  Wannier  functions  of  the  perfect  lat- 
tice. Furthermore,  as  the  localized  site  index  of  the  GWF's  moves  away  from 
the  defect  site,  the  GWF's  exponentially  approach  the  perfect  lattice  Wannier 
functions.  Then,  only  a few  perturbed  Wannier  functions  will  differ  from 
the  perfect  crystal  Wannier  functions.  Kohn  and  Onffroy  introduce  an  ab 
initio  method  for  calculation  of  the  GWF's  and  comnent  that  this  procedure 
is  easier  than  attempting  to  calculate  eigenfunctions  of  a non-periodic  sys- 
tem. This  theory  has  also  been  applied  to  a localized  impurity  state  between 
two  bands  and  the  determination  of  which  band  this  state  belongs  to."'1  Gay 
and  Smith"’'*''  have  extended  Kohn  and  Onffroy's  work  to  the  study  of  localized 
bulk  defects  in  three-dimensional  solids.  They  have  pointed  out  the  particu- 
lar usefulness  of  the  GWF  method  to  study  the  well  localized  deep  traps, 
complementing  the  effective  mass  formalism  for  shallow  impurity  levels  in  a 
band  gap  material.  Lattice  vacancies  are  also  considered  in  terms  of  a de- 
fect model.  In  contrast  to  the  Koster-Slater  method  of  handling  localized 
impurity  problems,  the  GWF  method  does  not  require  the  knowledge  of  the  im- 
purity potential  and  the  full  set  of  energy  eigenvalues  and  eigenstates  of 
the  perfect  crystal.  Gay  and  Smith  also  show  the  ease  of  self-consistently 
calculating  the  charge  densities  and  local  density  of  states  (LDOS)  knowing 
only  the  GWF's.  Since  the  charge  density  changes  in  the  region  around  a 
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defect,  self-consistency  has  been  shown  to  be  of  critical  importance.*11'7 
In  terms  of  the  bond  concept  in  solids,  the  GWF's,  due  to  their  orthonorma- 
lity, decay  smooth, y in  regions  where  no  other  functions  would  exist.  At 
a surface,  the  decaying  tails  of  the  GWF's  form  the  dangling  bonds  of  sur- 
face chemistry. 

Studies  of  the  application  of  GWF's  to  crystal  surfaces  have  been  done 
by  Rehr  and  Kohn248  and  by  Smith  and  Gay.24<3«250  The  importance  here  is 
the  detailed  characterization  of  the  electronic  structure  of  the  surface  in 
the  chemisorption  process.  GWF's  are  applicable  to  surfaces  because  of  the 
ample  experimental  evidence  that  the  strong  variation  of  electrical  proper- 
ties in  the  surface  region  decays  rapidly  into  the  bulk.  Describing  the 
surface  electronic  structure  in  terms  of  orthonormal  local  functions  which 
can  be  calculated  directly  (GWF's)  is  surely  an  advantage  over  using  the 
delocalized  functions  for  localized  surface  states. 

Smith  and  Gay249*280  test  the  GWF  formalism  on  exactly  soluble  one- 
dimensional Mathieu  potential  surface  problems.281"283  The  type  analytic 
model  potential  that  was  used  is  shown  in  Figure  34.  The  GWF  trial  functions 


Figure  34.  Model  Surface  Potential  Plot.  The  inset  shows 
the  relevant  part  of  the  energy  spectrum  of  the 
potential,  (after  Gay  and  Smith280’284) 
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were  chosen  as  modified  gaussiaris  with  a variational  parameter  and  are  shown 
for  the  first  three  lattice  sites  in  Figure  35.  This  shows  dramatically 
how  the  GWF's  rapidly  approach  the  symmetric  Wannier  functions  of  the  bulk 
and  the  smooth  exponential  decay  of  the  surface  site  GWF  into  the  vacuum. 


Figure  35.  Generalized  Wannier  Functions  Computed  with 
Gaussians  for  the  First  Three  Lattice  Sites. 
The  rapid  approach  to  the  symmetric  hulk 
Wannier  function  by  the  third  site  is 
illustrated,  (after  Gay  and  Smith750) 


The  comparison  of  the  exact  wave  functions  to  the  approximate  wave  functions 
calculated  from  the  GWF's  is  shown  in  Figure  36,  illustrating  the  accuracy 
of  results  using  simple  trial  functions.  The  theory  for  three-dimensional 
surfaces  and  their  complications  is  discussed  in  reference  246.  The  use 


96 


Figure  36.  Comparison  of  Approximate  with  Exact  Wavefunctions 
for  the  Mathieu  Surface  Potential.  The  approxi- 
mate wavefunctions  are  calculated  from  GWF's.  The 
top  three  plots  compare  wavefunctions  near  the 
(a)  top,  (b)  middle,  and  (c)  bottom  of  the  lowest 
band.  The  bottom  plot  compares  the  approximate 
and  exact  surface  state  wavefunctions.  (after  Gay 
and  Smith?50) 
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Needed  work  has  not  been  done  on  the  practical  application  of  General- 
ized Wannier  Functions  applied  to  two-  and  three-dimensional  models,  es- 
pecially involving  surfaces  and  interfaces.  Some  work  has  been  done  by 
Pei  and  Kohn  for  calculating  interface  energies  concerning  stacking  faults 
in  nickel,  resulting  in  practical  applications  to  transition  metals. ?<3r) 

In  a private  communication  with  J.  R.  Smith,  he  suggested  the  feasibi- 
lity of  applying  the  method  to  various  semiconductor  surfaces,  in  particular, 
silicon  with  its  bulk  diamond  structure.  Calculations  for  various  gallium 
arsenide  surfaces  also  seem  feasible.  Experimental  evidence  by  photoemis- 
sion and  low-energy  loss  spectroscopy  has  generally  shown  the  existence 
of  surface  states  localized  predominantly  on  the  arsenic  atoms.  Calculating 
the  GWF's  for  the  polar  faces  (either  arsenic-rich  or  gallium-rich  surface) 
and  the  subsequent  calculation  of  the  charge  density  and  LDOS  at  surface 
sites  would  help  corroborate  experimental  results  and  contribute  to  theore- 
tical predictions. 

The  effort  that  has  been  undertaken  in  the  study  of  the  Wannier  function 
techniques  has  shown  that,  until  recently,  the  technique  has  been  dormant 
for  several  years  --  overshadowed  by  extensive  band  structure  calculations. 
The  concepts  behind  the  theory  suggest  that  the  GWF's  are  a logical  exten- 
sion to  the  non-orthogonal  atomic  functions  while  preserving  the  idea  of 
Bloch  functions;  that  is,  explaining  the  nature  of  electrons  in  crystals. 

The  ability  to  calculate  GWF's  from  an  ah  initio  procedure  is  superior  to 
any  empirical  band  structure  procedure.  Also,  the  ability  to  calculate 
various  parameters,  even  the  Bloch  functions,  directly  from  Wannier  func- 
tions, reflects  its  generality.  Wannier  functions  normally  have  been  most 
useful  for  wide  band  gap  materials,  such  as  insulators,  but  have  been 
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applied  to  surfaces,  incorporating  strong  and  weak  perturbations  there. 

The  essential  ingredients  of  the  theory  show  that  it  can  effectively  be  ap- 
plied to  any  system  which  possesses  spatially  localized  states.  The  realm 
of  III-V  compound  semiconductors,  such  as  GaAs,  fits  this  application,  as 
well  as  various  interfacial  structures,  such  as  oxygen  chemisorption  on 
GaAs.  Photoemission  results  have  clearly  shown  the  retention  of  interface 
states  at  the  GaAs/oxygen  interface.  It  is  hoped  that  calculations  such 
as  the  ones  mentioned  here  will  provide  a theoretical  basis  for  practical 
interface  processes,  such  as  the  passivation  of  elemental  and  compound  semi- 
conductors. 

8.5.3  Localized  Electronic  Structure  Calculations  Using  Pseudo- 
potentials 

In  an  attempt  to  understand  the  behavior  of  electrons  in 
crystals,  the  pseudopotential  method  was  developed  as  a modification  of 
Herring's  orthogonal  ized  plane  wave  (OPW)  formal  ism. ?<5r’  The  method  relies 
on  the  assumption  that  the  effective  potential  for  the  core  electrons  may 
be  separated  from  that  of  the  valence  electrons.  This  is  in  line  with  the 
view  that  the  valence  electrons  are  solely  responsible  for  all  of  the  ord- 
inary chemical  and  physical  properties  of  the  system.  The  pseudopotential, 
then,  seeks  to  describe  the  interaction  of  a valence  electron  with  the  core 
ion  and  the  surrounding  cloud  of  other  valence  electrons.  Furthermore,  in 
contrast  to  the  view  of  the  Hartree-Fock  self-consistent  field,  a wide  range 
of  solid  state  properties  may  be  calculated  directly  using  the  pseudopoten- 
tial. Thus,  results  of  pseudopotential  calculations  may  readily  be  compared 
to  experimental  results  which  then,  in  turn,  could  be  used  to  modify  and 
refine  the  original  pseudopotentials,  i.e.  the  interaction  of  the  valence 
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electron.  The  pseudopotential  can  be  thought  of  as  inherently  containing 
the  effect  of  the  core  structure  while  retaining  the  detailed  structure  of 
the  valence  electrons.  The  lowest  states  of  the  pseudopotential,  therefore, 
correspond  to  the  lowest  states  of  the  valence  electrons.  The  pseudopoten- 
tial models  the  crystal  structure  as  a weak  periodic  perturbation  and  allows 
the  valence  structure  to  be  considered  by  the  well  known  nearly-free  elec- 
tron (NFE)  formalism.  Band  structures  are,  therefore,  easily  obtained. 

Rapid  theoretical  progress  coupled  with  applications  emerged  in  the 
early  1960's,  notably  by  Austin,  Heine,  and  Sham,?rj7  who  recasted  the  theory 
into  the  form  which  essentially  prevails  today.  This  era  ended  in  the  fa- 
mous paper  by  Cohen  and  Bergstresser258  in  which  they  documented  and  illust- 
rated pseudopotential  results  for  a total  of  fourteen  elemental  and  compound 
semiconductors.  Several  of  the  bulk  band  structures  were  found  to  be  well 
correlated  to  optical  reflectivity  and  photoemission  data  existing  at  that 
time.259 

Harrison's  book  in  1966  reviewed  and  contributed  to  the  pseudopotential 
theory,  particularly  for  metals,  up  to  that  time.280  It  was  essentially  re- 
placed by  the  authoritati ve  reference  by  Heine,  Cohen,  and  Weaire  in  1970. 281 
Most  present  day  work  is  based  on  the  foundations  laid  in  this  book  which 
includes  reviews  of  the  entire  theory  relating  to  the  various  types  of  pseudo- 
potentials in  use  as  well  as  the  detailed  fitting  of  pseudopotentials  to 
experimental  data.  Also  a large  section  on  cohesion  and  structure  is  in- 
cluded. 

It  is  worthwhile  to  emphasize  that  in  most  applications  of  the  pseudo- 
potential technique,  persisting  to  the  present  day,  the  pseudopotentials 
are  empirically  derived.  The  reason  for  this  is  that  considerable  success 
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has  been  achieved  in  fitting  pseudopotentials  to  one  experiment  and  then 
using  the  results  on  another  type  of  experiment.  Reference  261  cites 
several  instances  of  this  occurrence. 


In  more  recent  times,  a useful  self-consistent  pseudopotential  method 
for  handling  surfaces  has  developed. In  the  model,  periodicity  is  re- 
tained parallel  to  an  interface  and  electronic  screening  of  the  interface 
is  properly  taken  account  of  by  self-consistency.  The  partial  periodicity 
allows  one  to  apply  standard  pseudopotential  theory  as  outlined  in  refer- 
ence 261.  A particular  surface  structure  is  built  up  by  considering  a 
unit  cell  as  an  infinite  "slab"  of  the  crystal,  i.e.  a structure  possessing 
a surface  on  each  side  of  the  slab.  As  the  slab  is  periodically  repeated, 
the  surface  layers  are  built  maintaining  a two-dimensional  periodicity. 

One  need  only  to  consider  enough  layers  so  that  the  two  surfaces  of  the 
periodically  repeated  system  do  not  appear  to  affect  one  another,  i.e.  by 
allowing  enough  layers  for  one  surface's  effect  to  decay  into  the  "bulk". 

The  self-consistent  pseudopotential  method  described  here  requires  three 
inputs:  (1)  the  empirical  pseudopotential,  V^p,  usually  constructed  from 
fitting  bulk  energy  spectra;  (2)  the  geometrical  arrangement  of  the  slab; 
and  (3)  atomic  ion  core  pseudopotentials,  V^.  The  parameters  of  are 
fit  to  a Heine-Abarenkov  core  potential  as  outlined  in  reference  261. 

The  entire  procedure  will  now  be  outlined  briefly  by  resorting  to  the 
block  diagram  in  Figure  37.  VpMp  initiates  the  self-consistent  loop.  The 
total  charge  density,  p (r ) , is  then  calculated  followed  by  a Hartree  screen- 


ing potential,  V^,  and  a Slater  statistical  exchange  potential,  Vx>  and 


Vx  are  added,  resulting  in  the  screening  potential,  VSCR, 


which  is  then  added 


to  V jgN  to  yield  a total  pseudopotential,  V^.  VT  is  then  retained  for 
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comparison  to  the  next  iteration  until  self-consistency  is  achieved. ?r,? 


STEPS  IN  ACHIEVING  SELF-CONSISTENC Y 


Figure  37.  Block  Diagram  Indicating  the  Basic  Steps  in 
Obtaining  a Self-Consistent  Psendopotential 
(after  Schliiter,  et  al.?fi?). 


The  Si (111)  ideal,  relaxed  (inward  0.33  A),  and  (2x1)  reconstructed 
surfaces  were  the  first  applications  of  this  method  to  a semiconductor  sur- 
face. ?f,?  The  ideal  and  relaxed  surface  models  show  general  agreement 

with  theoretical  results. 2r,r>  We  illustrate  Cohen  and  his  group's  calcula- 
tion of  the  local  density  of  states  (LDOS)  and  total  charge  of  the  relaxed 
Si (111)  12-layer  slab  surface  model  in  Figure  38.  Layer  1 is  deep  within 
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the  bulk  and  is  typical  of  the  silicon  bulk  density  of  states.  Surface 
features  reveal  themselves  as  the  top  slab  (layer  6)  is  approached  showing 
the  prominent  Si (111)  dangling  bond  state  surface  state.  The  other  shaded 
peaks  are  correlated  to  the  four  back-bond  surface  states.  The  Si (2  x 1) 
reconstructed  surface  model  was  found  to  give  splitting  of  the  dangling 
bond  state  into  two  bands,  consistent  with  photoemission  data  and  other 
theoretical  results.265 


Figure  38.  Local  Density  of  States  for  a Relaxed  Si (111) 
Surface.  The  most  prominent  surface  states 
are  indicated.  The  total  valence  charge  as 
a function  of  the  coordinate  z perpendicular 
to  the  surface  is  shown  for  reference,  (after 
Louie,  et  al  .7F,t' ) . 
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The  electronic  structure  of  GaAs(llO)  ideal  and  relaxed  surfaces  has 
also  been  studied  using  the  pseudopotential  method. ?F,‘1  »2G7 »?68  For  the 
ideal  surface.  Figure  39  shows  the  total  valence  charge  density  terminated 
for  each  type  of  atom.  Figure  40  shows  the  charge  densities  for  the  dang- 
ling bond  surface  states  and  Figure  41  shows  the  LDOS  with  a corresponding 
table  correlating  the  positions  of  the  identified  surface  features  (shaded) 
with  respect  to  the  top  of  the  bulk  valence  band.?G7  The  fact  that  recent 
photoemission  data  does  not  agree  well  with  the  location  of  the  anion  and 
cation  derived  surface  states,  a relaxed  GaAs(llO)  surface  was  considered. 
Figure  4?  shows  the  LDOS  obtained  by  the  pseudopotential  calculations  for 
this  structure. 


Figure  39.  Total  Valence  Charge  Density  for  the  GaAs(llO)  Surface 
Plotted  in  the  (110)  Plane  Terminated  on  the  (a)  Ga 
Atom  and  (b)  As  Atom.  The  charge  density  has  been 
normalized  to  one  electron  per  unit  cell  volume 
^cen=8^2  A3.  The  contour  spacing  is  in  units  of  0.35. 

(after  Chelikowsky  and  Cohen767). 
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Go  Go  Go 


Figure  40.  Charge  Density  of  the  GaAs  (110)  Surface  for  the 
(a)  Ga  Da  rig  1 ing-Bond  and  (b)  As  Dangling-Bond 
Surface  States  in  the  Same  Plane  and  Normalization 
as  in  Figure  39.  The  contour  spacing  is  in  units 
of  2.0.  (after  Chelikowsky  and  Cohen?r,/)- 
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GoA»  (110)  SURFACE 
Local  daniity  of  itatet 


Feature 

Energy  (aV) 

Aa  5-like 

-9.0 

Oa  s-llke 

5. 70 

Aa  back  bond 

-2.0 

As  parallel  bond 

-0.8 

Aa  dangling  bond 

-0. 28 

Os  dangling  bond 

+ 1.0 

Figure  41.  Local  Density  of  States  for  the  GaAs  (110)  Surface 
(in  arbitrary  units).  The  total  density  of  states 
was  obtained  by  summing  over  all  layers.  Prominent 
surface  features  are  shaded  in  the  spectra  and 
identified  in  the  table,  (after  Chelikowsky,  et  al.?fi;)- 
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GaAi  (110)  Surface 


ENERGY  (.VI 


Figure  42.  Local  Density  of  States  for  a r 'Taxed  (110) 
Surface  of  GaAs.  The  LDOS  at  the  surface 
layer  for  the  ideal  case  is  also  displayed, 
(after  Louie,  et  al.?f>1*). 


The  effect  of  the  relaxation  is  shown  to  move  the  filled  As  surface  states 
down  about  1 eV  and  to  disperse  the  empty  Ga  surface  state.  These  results 
are  more  compatible  to  partial-yield  photoemission  data  than  the  results 
for  the  ideal  surface.770  Surface  states,  however,  still  are  found  in 
the  energy  gap  which  is  in  contradiction  to  experimental  evidence.  It  is, 
of  course,  suggested  that  a more  suitable  geometry  could  be  chosen  that 
would  improve  results. 

A pseudopotential  study  of  metal -semiconductor  interfaces  has  been  under- 
taken by  Cohen  and  his  group  using  the  jellium  model  to  describe  the  metal. 7r''* 
M/Si,  M/GaAs,  M/ZnSe,  and  M/ZnS  models  have  been  studied  because  of  their 
increasing  ionicity.  Pseudopotential  results  show  that  the  dangling  bond 
surface  states  are  removed  by  the  presence  of  the  metal,  but  are  replaced 
by  metal -induced  gap  states  (MIGS).  In  contrast  to  the  Bardeen  model,  in- 
trinsic surface  states  are  shown  not  to  play  a dominant  role  in  determining 
the  barrier  height  of  a Schottky  interface.  Calculated  barrier  heights, 
due  to  the  MIGS,  are  found  to  be  in  good  agreement  with  experimental  values. 

A simple  model  is  postulated  that  quantitatively  correlates  the  barrier 
heights  for  different  metals  on  semiconductors  of  varying  ionicity.771 

It  is  noteworthy  to  mention  that  recent  theoretical  work777  using  an 
empirical  tight-binding  Hamiltonian  for  GaAs(llO)  corroborates  with  most 
of  the  self-consistent  psendopotential  calculations.  The  pseudopotential 
method  has  also  been  recently  applied  to  GaAs-AlAs  heterostructures  [which 
can  be  grown  by  molecular-beam  epitaxy  (MBE ) ] with  new  fitting  parameters 
for  Ga , A1 , and  As. 773 

In  conclusion,  the  pseudopotential  method  has  recently  been  the  sub- 
ject of  several  investigations  on  surface/interface  electronic  states. 
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Calculations  have  shown  surface  features  that  are  in  relatively  good  agree- 
ment with  experimental  data  concerning  band  structure.  Also,  a model  has 
been  suggested,  based  on  pseudopotential  results,  which  puts  the  explana- 
tion of  Schottky  barriers  on  a firmer  footing.  As  experimental  data  be- 
comes more  accurate,  the  pseudopotentials  can  likewise  be  refined.  As 
results  of  pseudopotential  studies  advance,  insight  should  be  seen  into 
the  problem  of  passivation  of  semiconductors,  particularly  GaAs.  The  pre- 
sent requirements  of  two-dimensional  periodicity  and  using  identical  slabs 
may  turn  out  to  be  too  stringent  to  understand  the  details  of,  say,  oxygen 
chemisorption.  The  continued  analysis  of  intrinsic  surface  states  by  pseudo- 
potentials, using  better  geometrical  models,  may  presently  prove  to  be  the 
most  useful  effort  in  this  area. 
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I X . SURFACE  I NSTRUMENTATION : 

The  principal  analytical  techniques  used  durinq  the  program  were  Auger 

/ 

electron  spectroscopy  (AES)  and  X-ray  photoelectron  spectroscopy  (XPS).  The 

/ 

advantages  of  high  resolution  AES  and  XPS  were  pointed  out  earlier  in  this 
report.  Most  of  the  facilities  required  to  conduct  the  program  were  avail- 
able in  laboratories  of  the  Applied  Sciences  Laboratory  of  the  Engineering 
Experiment  Station  or  in  the  School  of  Physics.  The  high  resolution  AES 
and  XPS  studies  were  performed  during  visits  to  the  Applications  Laboratory 
of  Physical  Electronics  Industries  in  Eden  Prairie,  Minnesota. 

At  the  beginning  of  the  program  in  1975  the  ultra-high  vacuum  system 
that  was  available  provided  combined  LEED  and  Auger  capabilities  with  a re- 
tarding grid  analyzer,  sputter-ion  gun  and  a quadrupole  residual  gas  analy- 
zer. This  system  was  also  equipped  with  a turbo-molecular  pump,  a titanium 
sublimation  pump,  a large  ion  pump  and  a gas  manifold. 

During  the  first  year  system  modifications  were  made  to  improve  the 
detection  sensitivity  of  the  Auger  technique.  These  modifications  included 
the  addition  of  a PHI  Model  10-155  Cylindrical  Mirror  Analyzer,  a Model 
11-010  Electron  Gun  Control  Unit  and  a Model  11-500A  Auger  System  Control 
Unit. 

During  the  second  year  a PHI  Model  04-191  Sputter  Ion  Gun  and  a Model 
20-115  Control  Unit  were  added  to  provide  higher  ion  energies  and  a raster 
capability  for  more  uniform  sputter  etching  of  larger  areas.  Also  an  auxi- 
liary titanium  sublimation  pump,  adapted  from  an  Electro-Ion  pump,  was 
installed  directly  on  the  chamber  for  increased  pumping  speed.  Improved 
data  handling  was  accomplished  by  using  a Nicolet  Model  535  Signal  Averager, 
purchased  under  the  contract.  This  unit  has  been  particularly  useful  for 
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increasing  the  signal-to-noise  ratio  of  weak  Auger  derivative  signals  without 
sacrificing  resolution  and  it  has  been  used  to  obtain  N(E)  data  for  energy 
loss  measurements.  Fabrication  of  MOS  and  Schottky  barrier  test  structures 
within  the  Auger  system  was  made  possible  by  using  an  evaporator  mounted  on 
the  rear  flange  of  the  system.  This  in  situ  evaporator,  pictured  in  Figure 
14  of  our  Interim  Report  dated  23  January  1977,  provided  a capability  for 
aluminum  or  gold  evaporation  with  appropriate  shuttering  and  thickness  con- 
trol with  a quartz  crystal  thickness  monitor. 

During  the  third  year  another  evaporator  was  designed  for  :nclusion  in 
the  system  when  the  LEED  optics  was  also  in  place  on  the  back  flange.  With 
this  arrangement  surface  characterization  of  GaAs  samples  before  metal  depo- 
sition can  include  crystallographic  information  from  LEED  as  well  as  the 
measurement  of  surface  contamination  and  stoichiometry  from  AES. 

The  current  LEED/Auger  ultra-high  vacuum  system  is  shown  in  the  photo- 
graph of  Figure  43. 

The  majojjty  of  the  improvements  in  the  instrumentation  for  surface 
studies  can  be  attributed  directly  to  this  ONR-sponsored  program.  We  anti- 
cipate acquiring  in  the  near  future,  under  a National  Science  Foundation 
equipment  grant,  a multiple  technique  surface  analysis  instrument  which  will 
include  high  spatial  resolution  ( ^0.2  pun)  SAM,  a high  energy  resolution 
analyzer  for  AFS  and  XPS,  and  a SIMS  attachment.  This  equipment  will  greatly 
aid  investigations  of  the  type  described  in  this  report. 
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Figure  43.  Surface  Analysis  Instrumentation  Used  for  LEED, 
AES,  LELS  and  Other  Surface  Studies. 


X.  IDENTIFICATION  OF  PERSONNEL: 

The  following  personnel  participated  in  the  program  during  the  period 
1 May  1975  to  30  April  1978. 

E.  J.  Scheibner,  Research  Professor  of  Physics  and  Head  of  the  Surface 
Physics  and  Chemistry  Laboratories  of  the  Solid  State  Sciences  Division,  was 
Principal  Investigator  of  the  proqram.  He  received  his  Ph.D.  in  Physics 
from  Illinois  Institute  of  Technology  in  1955,  has  an  M.S.  degree  in  Physics 
from  Illinois  Tech  and  a B.S.  deqree  in  Physics  from  Georgia  Tech.  Dr. 
Scheibner  is  primarily  responsible  for  surface  analytical  instrumentation 
including  low  energy  electron  diffraction  (LEED)  and  Auger  electron  spec- 
troscopy (AES).  He  has  extensive  experience  in  thin  films,  semiconductors, 
ultra-high  vacuum  and  surfaces.  His  current  research  interests  relate  to 
surface  or  interface  effects  in  solid  state  devices  and  applications  of  AES. 

He  regularly  teaches  a graduate  course  on  surface  physics  as  part  of  an 
interdepartmental  academic  program  in  surface  science  and  technology,  and 
directs  the  thesis  research  of  Ph.D.  students  in  Physics. 

M,  W.  Ribarsky,  a Research  Scientist  in  the  School  of  Physics,  was  in- 
volved in  the  theoretical  phases  of  the  program,  particularly  the  analysis 
of  Auger  and  X-ray  photoelectron  spectra  and  the  development  of  methods  for 
calculating  energy  levels  for  clusters  of  atoms  at  the  surface.  Dr.  Ribarsky 
received  his  Ph.D.  in  Physics  from  the  University  of  Cincinnati  in  1974  and 
his  M.S.  and  B.S.  in  Physics  from  the  same  institution  in  1970  and  1967, 
respectively. 

Wal ter  H.  Hick! in,  Research  Engineer,  had  direct  responsibility  for 
operating  and  maintaining  the  LEED/Auger  equipment.  He  is  highly  competent 
in  the  operation  of  several  other  analytical  instruments  including  the  scanning 
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electron  microscope,  transmission  electron  microscope,  and  X-ray  and  electron 
diffraction  units.  In  addition,  he  has  experience  in  producing  films  of  a 
variety  of  metals,  insulators  and  semiconductors  by  evaporation,  sputtering, 
electroplating,  chemical  vapor  deposition  and  molecular  beam  epitaxy  tech- 
niques. 

James  M.  Hill,  III,  a graduate  student  in  Physics,  was  involved  in  the 
early  work  on  GaAs  surfaces.  Mr.  Hill  received  his  M.S.  Degree  in  Physics 
in  December  1976  and  accepted  employment  in  Texas. 

Richard  Marucci,  a student  in  Electrical  Engineering,  was  involved  in 
measurements  on  silicon  MOS  devices  and  in  the  development  of  Auger  instru- 
mentation. Mr.  Marucci  received  his  B.S.  degree  in  Electrical  Engineering 
in  December  1976  and  has  accepted  an  electronics  position  in  another  division 
of  the  Engineering  Experiment  Station. 

Charles  S.  Brown,  a graduate  student  in  Physics,  was  involved  in  the 
preparation  and  analysis  of  thin  oxide  films  on  tantalum,  aluminum  and  gal- 
lium arsenide  by  optical  methods.  Mr.  Brown  successfully  passed  the  compre- 
hensive examination  in  the  Spring  of  1977  and  is  now  engaged  in  his  Ph.D. 
thesis  research. 

A.  Craig  Kenton,  a graduate  student  in  Physics,  was  involved  in  Auger 
studies  of  GaAs  surfaces  and  in  the  analysis  of  Schottky  barriers  on  GaAs. 

He  also  initiated  theoretical  studies  of  surface  states  on  GaAs  by  cluster 
methods.  Mr.  Kenton  successfully  passed  the  comprehensive  examination  in 
the  Spring  of  1977  and  is  now  engaged  in  his  Ph.D.  thesis  research. 


I 
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X I . PRESENTATIONS,  REPORTS  AND  PAPERS: 

E.  J.  Scheibner,  "Solid  State  Reactions  on  Semiconductor  Surfaces," 

presented  at  Workshop  on  Solid-Phase  Epitaxial  Reactions  held  in 
Hughes  Research  Laboratories,  Malibu,  California,  August  8,  1975. 

E.  J.  Scheibner,  "Studies  of  Solid  State  Reactions  on  Elemental  and 
Compound  Semiconductor  Surfaces,"  Interim  Technical  Report, 

Contract  No.  N00014-75-C-0909  with  Office  of  Naval  Research, 

January  21 , 1 976. 

E.  J.  Scheibner,  "Interface  Structures  in  MOS  Devices  and  Their  Effect 
on  Reliability,"  presented  at  the  meeting  of  the  International 
Society  of  Hybrid  Microelectronics,  Huntsville,  Ala.,  January  28, 
1976. 

M.  W.  Ribarsky  and  E.  J.  Scheibner,  "Analysis  of  High  Resolution  Auger 

and  X-ray  Photoelectron  Spectra  from  Thin  Oxide  Films  on  Aluminum," 
presented  at  Amer.  Phys.  Soc.  meeting  in  Atlanta,  March  20  - April 
1,  1976,  Bull.  Amer.  Phys.  Soc.,  21,  432  (1976). 

E.  J.  Scheibner  and  W.  H.  Hicklin,  "Characterization  of  Metal-Oxide 
Systems  by  High  Resolution  Electron  Spectroscopy,"  presented  at 
the  30th  Annual  Symposium  on  Frequency  Control , Atlantic  City, 

N.J.,  2-4  June  1976;  Proceedings  30th  Annual  Symposium  on  Frequency 
Control,  U.  S.  Army  Electronics  Command,  Fort  Monmouth,  N.J.  (1 976; 
pp.  240-248. 

E.  J.  Scheibner,  "Studies  of  Solid  State  Reactions  on  Elemental  and 

Compound  Semiconductor  Surfaces,"  Interim  Technical  Report,  Con- 
tract No.  N00014-75-C-0909  with  Office  of  Naval  Research,  January 
28,  1977. 

Y.  E.  Strausser,  E.  J.  Scheibner  and  J.  S.  Johannessen,  "Observations 
of  A1 2 O3  and  Free  Si  in  the  Interface  Between  A1  Films  and  Si02," 
accepted  for  publication  in  Thin  Sol  id  Films. 

M.  W.  Ribarsky  and  E.  J.  Scheibner,  "Analysis  of  Auger  and  X-ray  Photo- 
electron Spectra  with  Application  to  Thin  Oxide  Films  on  Aluminum," 
unpublished  report. 
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XII.  SUMMARY  AMD  RECOMMENDATIONS: 


The  physics  and  chemistry  of  interfaces  in  composite  planar  structures 
on  silicon  and  gallium  arsenide  single  crystals  were  investigated.  In  the 
study  of  interface  chemistry  the  direction  of  oxidation-reduction  reactions 
was  predicted  from  the  bulk  thermodynamic  values  of  the  free  energy  chanaes 
on  forming  the  various  inorganic  oxides.  Thus,  observations  were  made  of 
the  reduction  of  Si 0^  on  silicon  by  A1  to  form  Al^O^  and  the  reduction  of 
Ga^O^  by  A1  to  form  Al^O^.  Other  metals  were  predicted  to  behave  similarly 
to  Al.  The  identification  of  the  separate  solid  state  phases  in  the  compo- 
site structures  was  made  by  high  energy  resolution  Auger  electron  spectros- 
copy (AES)  and  X-ray  photoelectron  spectroscopy  (XPS)  techniques.  Further 
studies  of  this  type  are  recommended.  In  particular,  investigations  should 
be  made  of  oxidation-reduction  reactions  involving  metals  other  than  Al . The 
deoendence  of  the  width  of  the  reaction  zone  on  initial  oxide  thickness  and 
the  elemental  distribution,  e.q.  Si,  Ga  or  As,  within  the  metal  and  the  oxide 
should  be  determined. 

Experimental  MOS  and  Schottky  barrier  test  devices  were  fabricated  from 
planar  composite  structures  prepared  in  the  ultra-high  vacuum  system  and 
having  wel 1 -characterized  interfaces.  Because  of  the  importance  of  these 
type  devices  in  advanced  electronic  systems  further  efforts  should  extend 
the  fundamental  studies  towards  new  device  technologies  including  a considera- 
tion of  device  reliability  from  the  materials,  rather  than  the  statistical, 
viewpoint. 

The  attenuation  length  of  '1400  eV  electrons  in  thin  oxide  films  on 
aluminum  was  studied  with  high  energy  resolution  AES  and  XPS.  Theoretical 
interpretation  of  the  data  involved  the  removal  of  background  and  instrument 
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resolution  effects  and  separation  of  the  oxide  and  metal  peaks  into  Lorent- 
zian  and  Gaussian  components.  Measurements  of  the  integrated  intensities 
of  the  peaks  as  a function  of  oxide  thickness  fitted  to  predicted  curves 
gave  values  for  the  attenuation  length.  Since  relatively  little  data  exist 
for  electron  attenuation  in  oxides  it  is  strongly  recornnended  that  these 
studies  be  continued.  In  particular  the  importance  of  establishing  accurate 
sputtering  rates  for  oxides  should  be  noted  and  measurement  schemes  should 
be  devised  for  obtaining  attenuation  length  data  for  a range  of  electron 
energies. 

The  passivation  of  compound  semiconductor  surfaces  on  devices  has  proved 
to  be  much  more  difficult  than  that  of  silicon  surfaces,  yet  it  is  necessary 
for  achieving  stable,  reliable  devices.  We  have  reviewed  the  voluminous 
literature  on  the  subject  of  GaAs  surface  passivation  and  have  suggested  a 
possible  method  based  upon  the  solid  state  reactions  occurring  at  interfaces. 
The  main  point  is  that  attention  to  the  fundamental  physics  and  chemistry 
of  interfaces  is  likely  to  be  equally,  if  not  more  fruitful  than  the  previous 
empirical  approaches  to  passivation.  In  future  efforts,  consideration  of 
the  electronic  factors,  e.g.  the  distribution  of  interface  states,  which 
are  required  for  a good  passivation  method  should  be  stressed.  In  addition 
more  directed  studies  of  interface  chemical  reactions  should  be  carried  out. 

Localized  theoretical  models  for  predicting  the  electronic  structure 
of  groups  of  atoms  in  the  vicinity  of  interfaces  have  been  examined  briefly. 
The  particular  techniques  included  are:  (1)  molecular  cluster  calculations 
based  upon  an  unrestricted  Hartree-Fock  method;  (2)  the  pseudopotential 
method;  and  (3)  the  generalized  Wannier  function  method.  Data  resulting 
from  photoemission  and  energy  loss  spectroscopy  will  provide  the  experimental 
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verification  of  the  assumed  parameters  in  the  models.  For  example,  it  is 
expected  that  precise  atomic  geometrical  arrangements  will  be  required  in 
most  of  the  theoretical  calculations  to  give  good  agreement  with  the  distri- 
butions of  surface/ i nterface  states  obtained  experimentally.  However,  a 
sensitivity  analysis  will  have  to  be  made  for  each  of  the  theoretical 
methods.  It  is  strongly  reconrnended  that  studies  of  the  localized  electro- 
nic structures  be  continued  with  emphasis  cn  the  atomic  composition  and 
arrangement  at  interfaces. 
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APPENDIX  A 


GaAs  WET  CHEMICAL  CLEANING  PROCEDURE 

1.  Deionized  water  rinse 

2.  Blow  dry  - 30  P.S.I.  N? 

3.  Methanol  rinse 

4.  Vapor  degrease  in  trichloroethylene  - 2 minutes 

5.  Methanol  rinse 

6.  Deionized  water  rinse 

7.  Blow  dry  - 30  P.S.I.  N? 

8.  "Piranha"  bath  (1  part  deionized  water)  - 60  seconds 

(1  part  H?0?) 

(5  part  HjSOiJ 

9.  Hot  deionized  water  immersion 

10.  Deionized  water  rinse  - 1 minute 

11.  Blow  dry  - 30  P.S.I.  N? 

12.  Bromine  etch  (1  part  bromine)  - 60  seconds 

(75  parts  methanol ) 

13.  Methanol  rinse  - 60  seconds 

14.  Deionized  water  rinse 

15.  Blow  dry  - 30  P.S.I.  N^ 
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